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Aspects of the Turbulence Problem 
2 Survey Report 


By Hans W. Liepmann, Pasadena, Calif.2) 


« First Part 


I. INTRODUCTION 


y. It is an empirical fact that fluid flow at sufficiently high Reynolds numbers 
is always turbulent. For high Reynolds numbers the velocity components u, 
always depend explicitly on the time ¢, even if the boundary conditions are 
independent of time. No steady flow is possible. Furthermore, u,(x, t) is a 
- stochastic function in space and time. These properties define turbulent flow. 
Turbulence is a very general phenomenon. It exists in any type of fluid 
motion and appears to depend only upon a few very broad characteristics such as: 
(a) The existence of nonlinear or rather quasi-linear terms in the equations of 
9 - motion which represent the transport acceleration in Eulerian coordinates. 
_(b) The existence of a viscous damping term of an order higher than the non- 
____ linear terms. 
_(c) The three-dimensionality of the motion. 

The vast majority of the motions of fluids are turbulent. Consequently 
turbulence appears in many physical and astrophysical situations, sometimes 
as the prime mover and sometimes as a source of perturbation and of ‘“‘noise’’ 
“in a general sense. 

_ The theory of turbulence has the aim of determining the statistical proper- 
ties of u;(X, #), that is, the mean values of u,, the probability distributions, etc. 
A typical problem is, given the boundary or initial conditions of a flow such as 
the flow from infinity past a solid body, the flow through a pipe, etc., to find 
“the average forces, the average velocity field and its fluctuations. In general, a 
problem like this is still far from being solved. For example, it is impossible 
today to predict the flow in a simple diverging channel without recourse to 
empirical results. It is known empirically that the flow character changes 
‘radically if a certain limiting angle of divergence is surpassed. Turbulence 


theory should predict this angle. 
; 
3 1) California Institute of Technology. 
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The three most conspicuous steps in turbulence research in recent years are: ~ 
(1) The introduction of the concept of isotropic and homogeneous turbulence 

together with the development of a suitable formalism of description of — 

turbulence (G. I. TAyLor and von KARMAN). 

(2) The idea of local isotropy and the asymptotic behavior of turbulence at 
very high Reynolds numbers (KOLMOGOROFF, ONSAGER, and VON WEIZ- 
SACKER). 

(3) The discovery of the intermittent nature of turbulent shear flows (CORRSIN 
and TOWNSEND). 
(1) defined the simplest problem in turbulence for both theoretical and 

experimental research. Most of the recent ideas in turbulence have developed 

from the study of isotropic turbulence. 

(2) is a successful attempt to obtain asymptotically valid results for the 
behavior of the elements of turbulence of smallest dimensions which have 
resulted from a cascade process from larger elements. 

(3) is the discovery of the fact that in turbulent shear flow a type of super- 
structure appears which influences the transport properties of turbulence deci- 
sively. The importance of intermittency has not yet been fully exploited. It is 
probably one of the most important features of turbulent transfer problems and 
very likely holds the key to turbulent separation. 

Besides the central problem of computing the structure of the velocity field 
there are at least two more general problems in turbulence research: the 
question of the origin of turbulence and of transition from laminar to turbulent 
motion is one. The other is the effect of turbulent fluctuations on the motion 
of rigid bodies or on the propagation of waves. 

In the following no attempt is made to give a complete review of the tur- 
bulence problem. Indeed, any such attempt would require a monograph. 
Instead, an attempt is made to discuss some general aspects of the turbulence 
problem, in particular those which are often omitted in reviews of turbulence. 
The important results of isotropic turbulence research are only briefly and 
incompletely summarized. However, for isotropic turbulence there exists an 
excellent recent review of AGostin1I and Bass [1]!) as well as a number of 
earlier reviews [4], [6]. Stability theories and transition were omitted from the 
present review. The laminar stability theory can be found, for example, in the 
recent book of SCHLICHTING [11] and the review articles of PILLow and DryDEN 
[9], [5]. Thermal instability problems, which are most important for meteoro- 
logical and astrophysical applications, have been discussed by CHANDRA- 
SEKHAR [41]. 

Only the general characteristic features of turbulent shear flow are stressed. 
Consistent results on the detailed nature of turbulent shear flows have been 


1) Numbers in brackets refer to the Bibliography, page 341. 
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obtained only recently. To summarize these recent results and interpretations 
in a critical manner is a difficult task. Besides the space limitations, the author 
has not felt competent to undertake it here. For the details of shear flow 
reference is made to recent work of TOWNSEND [63], [65], LAUFER [52], Rotta 
[59], and KLEBANOFF and DIEHL [50]. 

In Section II a brief survey is given of elementary results of the theory of 
stationary, stochastic processes. This section is intended to introduce concepts 
like correlation function and power spectrum. It also provides a link between 
the terminology and methods in turbulence research and in the theory of 
stochastic processes. In Section III the response of linear systems to random 
forcing functions is discussed, and a number of examples are given where 
turbulence plays the role of the forcing function. This section is intended to 
indicate the variety of physical situations in which turbulence enters. The 
examples demonstrate that it is frequently possible to say something about the 
effect of turbulent fluctuations without a very detailed knowledge of the 
mechanism of turbulence itself. In Section IV the general problem of turbulent 
transfer is stressed and the concepts of local isotropy and of intermittency are 
discussed. F inally, a short review of isotropic, homogeneous turbulence is given 
in Section V. 

The author wishes to express his appreciation for critical discussions with 
Doctors J.D. Corr, P. A. Lacerstrom, and A. RosHKxo which played an 
important part in preparing this review. 


Dies A LON Ake LO CHS V1 CaP iO CESS ES 


As a preliminary in dealing with turbulence it is very helpful to discuss 
briefly some elementary results of stationary stochastic processes. For a detailed 
study of the subject the reader is referred to WIENER [18] and Rice [16]. 


1. Stochastic Function 


If J(t) is a stochastic function of ¢, usually the time, then J for any given ¢ 
is a random variable. The values of J at ¢,, ...,¢, form the components of a 
random vector. The probability structure of /(f) is determined completely by 
the joint probability for each such random vector. 


2. Mean Values 


The mean values of functions of J(¢) are computed from the probability 
distribution and hence are ensemble averages. For certain stationary processes 
an ergodic property exists and ensemble averages can be replaced by time 
averages. Averaging will be indicated by a bar over the functions, 
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3. Stationarity 


A stochastic process is called stationary if all mean values and probability 
distributions are independent of the choice of the origin of ¢. The process is thus 
homogeneous in time. 


4. Correlation Function, Power Spectrum 


J] (t) = constant by stationarity. It is convenient to normalize and put 


J) =0. (II-1) 
The correlation function or auto-correlation function y(t) of J(¢) is defined by 
o(t) = JO IE+2) = 9-2). (II-2) 


In the following it will always be assumed that y(t) + 0 for large 7 sufficiently 
rapidly to ensure the existence of the integrals defining the power spectrum 
@(w) and the time scale 0. 


ie.°) 


P(w) = a / y(t) cosw t dt, (II-3) 
y(t) = [ @P(w) coswmt da, (II-4) 
(0) 6 = [9 dc = = (0). (II-5) 


0 


(II-3) and (II-4) are often called the Wiener-Khintchine relations. In the 
theory of isotropic turbulence these relations were introduced independently 
by G. I. Taytor. 

The assumptions made exclude periodicities in J(t) and a discrete part in 


the power spectrum. For the study of turbulence this appears to be a natural 
assumption. 


5. Ergodic Property 


A stochastic process is called ergodic if the ensemble averages may be com- 
puted as time averages. A sufficient condition is the existence of a time scale 0 
and suitable similar integrability conditions for correlations of higher order. 
For example, consider the time average 


Tye rf{I® dt. (11-6) 
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I(T) is a stochastic function of T. Its mean value is equal to the mean value 
J (t). If it can be shown that the dispersion of /(T) becomes negligible as T 
becomes large then one measurement of /(7T) determines its mean and hence 
also the mean of J(¢). Since J(t) = 0 the dispersion of I(T) is given by 


| [TATE  dipdig = a ice dt, dt, (11-7) 


or 
r= & frames <d five (IL-8) 


Hence I(T)? +0 as Too. Similarily, to compute the auto-correlation 
function of /(¢) using the time average, an integrability condition for the 
quadruple correlation is needed, and so forth. In the special and important case 
where the mean square of a normally distributed /(¢) is computed from the 
time average the dispersion formula becomes quite similar to (II-8). Namely, if 


Dy, 
L(T) => | J2(t) dé (I1-9) 
then : 
[L(T) — L(T)}* = se / (T — 2) p(x) de. (II-10) 


0 


The intuitive idea underlying the above reasoning is simply the following: if a 
time 7, exists such that J(¢+ 7) is essentially independent of J(t) for + eda) 
then averaging over a time T > Ty, is equivalent to taking the ensemble 
average over 1 systems with 

Ui tod Sie . 
It is interesting to compare this reasoning with the proofs for the law of large 


numbers (for example [14]). 
6. Mean Values of the Derivatives of /(t). Spectral Moments 
Differentiating (II-2) with respect to t yields 


g@=JOTEHD=SE-DIO, | 
Seo Mmerie-staire | 


Hence 


(OR ie a i2 (11-12 
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and using (II-4) it follows that 


"(0) = — fe w* @(w) do. (11-13) 
0 


The derivatives of the correlation function at t= 0 are thus simply related 
to the moments of the power spectrum. In general we have 


(J)? = | co®™ Dl) dew = (— 12-4 p20), (11-14) 


0 


The existence of mean square derivatives of J thus implies the existence of the 
corresponding moments of the spectrum. If J has all mean derivatives then 
@(w) has to vanish rapidly enough for large w to guarantee that all moments 
exist. This result is applied later on in arguments concerning the spectrum of 
turbulence. 

Furthermore, it follows from (II-11) and (II-2) that 


Jj’ =0. (II-15) 
J(t) and its derivative at the same time have zero correlation. The same is 


obviously true for any correlation between derivatives differing by one order, 
for example 


J J’ =0, ete. (II-16) 


7. Probability Distribution 


To define J(t) it is not sufficient to give y(t) or B(w) ; it is also necessary to 
prescribe the probability distribution of J and the joint distributions of J(4), 
J (tz), etc. It is here that an essential difference between turbulence and linear 
random noise lies. In the theory of random noise, Brownian motion and related 
topics it is often possible to assume that J (t) is normally distributed (see, for 
example, RICE [16] and UHLENBECK and WANG [17]). 

If J and J’ have a Gaussian distribution then it follows from (II-15) that 


they are also statistically independent. The joint probability is the product 
of the distribution of J and J’ 


bold) = (2. p(0)]-"* exp [— 525], (I1-17a) 


A(F) = [—22 p"(0)]-"” exp| + yn | (II-17b) 


and 


Pall, J’) = bolJ) Pi(J’) - (II-18) 
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For more general cases (II-15) does not imply (II-18). If /(¢) denotes one 
component of the velocity fluctuations in turbulent flow, for example, it 
follows that because of the nonlinear terms in the equations of motion 


(7)? 4-0. (11-19) 


Consequently J’ can not have a Gaussian distribution in this case. 
The application of the central limit theorem to the probability distribution 
of /(t) is best seen by representing J (¢) as a Fourier-Stieltje’s integral 


+ ©Oo 


ee / ei! dF (2), (11-20) 


— oo 


where F(A) is a stochastic function of A. For stationarity it is necessary that the 
increments dF(A) be uncorrelated. This is seen by forming the correlation 
function 


ple) = JO TE+ a =f [ei Mte te dF (a) dF*(u). (01-21) 


—0o 


For a stationary process thus 


dF (A) dF*(u) = \dF (A) |? (2 — po) 


and 
pte) = fe dF). (11-22) 
Comparing (II-22) with (II-4) shows that 
@(A) da = \dF (A) |?. (II-23) 


If in addition to being uncorrelated the increments dF'(/) are independent one 
can apply the central limit theorem to (II-20), considering J(¢) for any ¢ as 
the limit of a sum of independent increments ~ dF (A). One concludes that J (t) 
is normally distributed. For the application of the central limit theorem it is 
essential that |dF(A)|? ~ da. 

For linear systems it is thus often possible to assume a normal distribution. 
This is done, for example, in problems of the shot effect, Brownian motion, and 
black body radiation. In homogeneous turbulence the increments dF(A) are 
uncorrelated, but not independent. The central limit theorem does not apply. 


8. Zeros and Extrema of a Stochastic Function 


Rick [16] has obtained a formula for the number of times per unit time 
N,(&) a stochastic function J] (t) passes through the value &. N,(é) can be ex- 
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pressed in terms of the joint probability 9, (&, 1) of the function J and its deriv- — 


ative J’ 
Nolé) = [ In| bulé, 0) dn. (11-24) 


(IJ-24) expresses the fact that J spends the time fo,(€, 1) d& dy between € and 
— + dé having a slope in the interval 7, 7 + dy. The crossing of the interval 


dé takes the time t = dé/ || and hence (II-24) results for the expected number ~ 


of crossings, regardless of slope. Applying (II-24) to the derivative of J an 
analogous result can be found for N,(7), etc. Specifically, No(0) gives the 
expected number of zeros of J(¢) per unit time, N,(0) the expected number of 
extrema. 

(II-24) is important when one is interested in the number of times J reaches 
or exceeds a certain value & rather than the length of time & is exceeded. 
Structural fatigue problems are typical problems of this type. 

For a double Gaussian distribution N,(0) is easily expressed as: 


72 11/2 —— 
No) == |4- |" = Vee. (11-25) 
Consider for the moment J as a velocity component, the independent variable ¢ 
as a space variable, and w as a wave number in space. Viscous dissipation of a 
wave in hydrodynamics is proportional to the square or mean square of the 
wave number. ([1-25) can then be interpreted in the following way: the mean 
viscous dissipation of a one-dimensional, Gaussian, stochastic velocity field is 
equal to the dissipation of a harmonic wave having the same average number 
of zeros. 

For a non-Gaussian process like turbulence (II-25) is not exactly correct. 
Measurements [91] show that the dissipation in turbulence is larger than the 
harmonic wave with an equal number of zeros. This is apparently due to the 
fact that the nonlinearity introduces a saw tooth character into the fluctuations. 
For the same number of zeros, the dissipation is then larger. 

Mathematically it means that fo,(&, 7) is not the product 4,(é) P1(7) and 
not double Gaussian. The step from (II-24) to (II-25) can not be made. 


9. Generalization to Homogeneous Fields in More Than One Dimension 


In the study of homogeneous turbulence, vectors or tensors that are sto- 
chastic functions of vectors are needed. The concept of power spectrum and 
correlation function is easily extended. For example, to study isotropic, homo- 
geneous turbulence, von KARMAN introduced the correlation tensors of various 
ranks. The correlation tensor g;;(7), for example, is defined in terms of the 


velocity components ; at two points P,(R) and P(R + 7) by 
isl?) = u(R) a(R +7). (II-26) 
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Generalizing ey and (II-4) we have 
®,,(R) = sae | dr gilt) e®", pli ey dK ®,,(K) e°'*", (11-27) 


There is no difficulty in extending this generalization further to tensors of 
higher rank and with more independent variables. In the following, for example, 
it will | be useful to consider correlation tensors in space and time. In this case 7 
and K are four-dimensional vectors. 

In general, all properties of stationary processes, that is, processes that are 
homogeneous in time can be generalized to random fields which are homogeneous 
in space or in space and time. In dealing with homogeneous fields we will thus be 
able to use an ergodic property under similar conditions to those in Section IT. 51). 
Evidently relations similar to (II-12), (II-13), and (II-14) hold by virtue of 
(II1-27) and suitable generalizations. 


fi URDULCENCE PREEECTS UPON LINEAR SYSTEMS 


In the following the response of linear systems to random forcing functions 
arising from turbulent fluctuations will be discussed. Turbulence acts as a 
random forcing function in a large variety of physical situations ranging from 
the motion of an airplane in gusty air to the scintillation of stars. Furthermore, 
the process of measuring turbulent fluctuations by means of a hot-wire anemo- 
meter or any other probe consists essentially in measuring the response of a 
system upon which the turbulent fluctuations act as an impressed force. 

Broadly speaking, turbulent fluctuations will have similar effects for large 
scale phenomena as molecular fluctuations have for small scale phenomena. 
Hence one encounters macroscopic phenomena similar to Brownian motion, 
Rayleigh type scattering for radio waves, etc. In the measuring process, 
however, the phenomena encountered in turbulence differ from similar obser- 
vations on molecular structures. For turbulence measurements it is possible 
in principle to keep the fluctuations within the measuring apparatus itself 
negligibly small. A separation between the measuring instrument and the 
fluctuations can be realized here. On the other hand, we are interested in much 
more detailed information concerning the fluctuations than is the case in most 
molecular fluctuation phenomena. 


1. Use of the Correlation Functions and Power Spectra 


We will assume in general that the fluctuations are homogeneous in time 
and space. Exceptions will be noted explicitly. 


1) Actually, it is possible in this case to compute ensemble averages by averaging over one 


space coordinate only. 


: 
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The correlation function appears in the discussion of the response of a 
system, due to the use of an integral representation and the averaging process. 
In one dimension let 


o(T) = if j(t) dt, (IIT-1) 
then 
»(T) = cf / jt) f(t) dt dt’. (III-2) 


Because of homogeneity 


f(t) Ht’) = p(t — #’) . 


where @ denotes the correlation function or ‘‘auto-correlation function” of f. _ 
Since depends only upon |t — ¢’| the integral can be reduced to a single in- | 


tegral 
; 
v2=2 / (T — ) y(n) dn. (III-3) 


“ 


0 


(III-3) has two limiting cases which are often useful. p(x) may be nearly equal 


to y(0) = f? in the interval (O, T). That is, the correlation extends for long 
times. Then 


v*(T) = f? T? (T<6)  (III-4) | 
or else p may drop to zero ina time T, < T. In this case put : 
p(n) = 0 d(n) (6 ~ T,) 
where 0 is the Dirac delta-function. Then 
vX(T) =272?TO. (TS 6)  (IIL-5) | 


Equations (III-4) and (III-5) are, of course, well known in diffusion theory if 
v?(T) is interpreted as the mean square displacement of a particle at the time 
T and f? its mean square velocity. In this case the velocity is homogeneous in 


time, the displacement not. Similar results can be written immediately in 
space or in space and time. 


2. Linear Langevin-Type Equations 


Consider now a linear operator {;; (for example, corresponding to a system 
of partial differential equations). Let v,(%) denote a vector field, /,(7) a forcing 
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function 
Lis(ds) = fi - (III-6) 


For a stochastic f,, (III-6) is a generalized Langevin equation. For an infinite 
domain 


v8) = | By) HE+9), (111-7) 


where y;;(s) denotes the fundamental solution of (11-6). For homogeneity it is 
necessary to assume constant coefficients in @;,. 

From (III-7) we obtain the mean square v,, or the correlation v,,(7) exactly 
as in the simplest case given above (III-2). Thus 


vi +1) vf) = ff a8 45" yal) Yom) AE +7 +5) P48). (ILS) 


In this way the correlation function of the response v,; can be expressed in 
terms of the correlation function of the forcing function f;,: 


vi(1) = a a dS 48’ Yes) Yim(S') fom” +S — Sy. (III-9) 


Introducing the power spectrum tensors, V;; and F,,, of v;; and f,;, and I’;; 
defined by?) 


Ty(R) = | a yi) (111-10) 
we have 
Vii(K) = Dn(K) I}, (K) Fin(B) « (11-11) 


7 and K do not need to be three-dimensional vectors. If 2;; includes time deriv- 


atives 7 will include time as one variable and K the time angular frequency w. 
If £,; does not involve time derivatives, time enters into the problem as an 
external parameter only. In this case (III-11) can be written with a three- 
dimensional k 


= sy => > 


Vis(b, t) = Dialh, 8) PAB, t) Fill, f) (111-12) 


(III-12) is similar to the equations given recently by UBERor and Kovasznay 
[33] in slightly different terminology and derivation. 

Equation (III-11) is the generalization of the relation linking the power 
spectrum ®(w) of a forcing function /(¢) to the power spectrum of a response 
v(t). In the latter case we have simply 


V(w) = |I'(@) |? ®(@) (III-13) 

1) The J7,; are used in ZEILON’s construction of the fundamental solution. See for example: 

FRANK and von Mises, Die Differential- und Integralgleichungen der Mechantk und Physik, Vol. 1. 
(F. Vieweg & Sohn, Braunschweig, 1930), p. 865. 


; 
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where 
or operator linking v(t) to f(é). Similarly J’;;(K K) enters into the relation between 
V,; and F,;, the corresponding power spectrum tensors. J’;, is obtained from 
the operator L;; in exactly the same way as the complex admittance is ob- 
tained from a one-dimensional operator: 0/0x; is replaced by « K; in &;;. 

In general, problems involving boundaries of the region are more compli- 
cated. As an example, let 2;; include time derivatives. Some interesting situa- 
tions lead to the problem where a forcing function /,(x, ?) which is homogeneous 
in space is switched on at time t = 0. The space correlation tensors v,,(7, t) are 
to be computed as functions of time. An example of this type will be discussed 
later on.) We have 


Qylv) =f 20); 659;50 SO), (III-14) 
Then 


t | 


0(%, 4) = [ ds [de yl — 8,t—0) f(s, 0) (III-15) 
e : : 


or 


=[ [a ds’ | 


Her . . Pera (III-16) 
Ke / / do do’ y,; (8, t — 0) Yim (s’, t — 0’) fim(y + s — s',o— 0"). | 
0 
The relation for the spectral tensors becomes 
V,,(k, t) ah / do do’ I'y(k, t — 6) D'jm(k, t — 0") Fim (R, 6 — 0’). (11-17) 


} 


Thus in this type of problem it is necessary to know the space-time correlation 
of the forcing function even if only the space correlation of the response is} 
desired. 


3. Some Typical Examples 


(A) Hot-Wire Anemometer 


Equation (III-6) describes the process of measurement of fluctuating phe- 
nomena in most instances. An elementary case is the response of a hot-wire 
anemometer exposed to the fluctuating velocity u(t). Here the voltage v(t) 
satisfies the approximate equation 


ai 
M“-+u=cult), (ITI-18) 
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where M denotes the so-called time constant of the hot-wire and c its sensitivity. 
Hence, the observation of the mean square voltage fluctuation yields 


2 @D ° 
Ppl i ee by aie i Dw) do. (III-19) 
0 


The so-called compensation method of hot-wire anemometry (DRYDEN and 
KUETHE [25]) consists in removing the factor 1+ M?w? from (III-19) by a 
second properly chosen linear circuit. Let v(t) denote the output of this second 
circuit 


QV) =. 
Then & has to be chosen such that |/(w) |? = 1 + M*? w®. Then 
Viwu?. 


For spectral measurements | J\(w) |* has to possess a very sharp maximum for, 
say, ® = W,. If a circuit possesses this quality, (III-13) gives for the mean 
square response 


0? = flo) | |T(@)|? deo. 


If the scale of the fluctuations is not very large compared to the wire 
length we have to take into account that the wire averages over its length J. 
Hence the wire does not measure u(¢) but rather 

1 


Ul) => 2 u(t, 2) dz. (ITI-20) 
0 

If only the mean square of w is required, or certain moments of the spectrum, 

the corrections are easily obtained. For example, we have in homogeneous 


turbulence (III-3) 
I 


at we [ (i — 2) o(2) dz. (III-21) 
0 
For short wires g(z) can be developed in powers of z and introducing the 
so-called microscale / of isotropic turbulence, we find 
pa — 1 [2 
aw lia a =h a], 
For a very long wire we may apply the equivalent of (III-5) and find 


A. 


ma 2A, 


where A now corresponds to the integral scale of turbulence. 


2 
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For the measurements of spectrum and correlation functions the length 
corrections are considerably more involved. UBEROoI and Kovasznay [33] have 
recently given a general approach to the problem using equations equivalent 
to (III-11) or (III-12). They found the correction for wire length in the measu- 
rement of the one-dimensional spectrum of isotropic turbulence: 


al yea) 2) 
Gly) == [ W (zVae— me) AOR) ae, (tL-22) 
ky : 
where G(f,) denotes the one-dimensional power spectrum as measured with a 
wire of length J. E(k) is the three-dimensional spectrum, k = VR + k3+ kz, 
and W a function defined by ; 


W(x) == (“S7\* (II1-23) 

In [33], W(x) is plotted and several special cases are discussed. An inter- 
esting result is obtained if the true spectrum is assumed to be of the form le 
In this case it is shown that the measured spectrum at high frequencies approa- 
ches the form k,"*”), This is an important result with respect to comparison 
of spectrum measurements with theoretical results concerning the form of the 
spectra where the exponent 7 is often of primary concern. 


(B) The Effect of Turbulence Upon the Propagation of Waves 


Turbulence in the atmosphere, in wind tunnels, etc. will cause fluctuations 
in the refractive index for the propagation of waves. The velocity fluctuations 
themselves as well as the temperature fluctuations will amount to variations 
in the propagation velocity for sound waves while the fluctuations in density | 
and temperature will cause fluctuations in the refractive index for electro-. 
magnetic radiation. Problems of this type have become rather important in 
recent years and a number of investigations have been carried out ({20] to [22], 
[27], [29] to [32]). Here a few typical effects will be discussed. We have to 
distinguish between two types of problems: the effect of the fluctuations may 
appear mainly as a “‘noise problem”’, that is, the interest may li in the primary 
signal while turbulence causes mainly a disturbance and diffusion of this signal. 
Or else the disturbance in a primary signal or wave may be used to obtain 
information about the turbulence pattern itself. Especially simple cases occur 
when either the wave length of the radiation is small compared to the scale 
of turbulence and geometrical optics applies, or the scale of turbulence is 
small compared to the wave length and we deal with a problem in diffraction. 
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The Passage of a Light Ray Through a Turbulent Boundary Layer 


As a simple problem consider the diffusion of a ray of light which traverses 
a boundary layer of thickness 6. Assume the incident ray normal to the layer, 
in the direction x. Within the layer the refractive index n is assumed to be 
of the form 


n= n(x) [1 + 14 (z, y, 2, 2)], 2 = n(x) . (III-24) 


The unit vector s in the direction of the light ray can be written 


onl, ea, 20), 


For small 1, and « it follows from the geometrical optical equation: 


as i 
g or a oid (IIT-25) 


A ni iig EH 


where 


= On, 
Vem =m — >, 4- 


(III-25) corresponds to diffusion in a field of force. The force field is due to 
the mean gradient of . Here 1, is homogeneous in space and time. a is homo- 
geneous in time, but not in space. Integrating, squaring and averaging of 
(I1I-25) yields 


6 
a? = [ng(d)]-* | | mol) molé) (V2 Male (V2 ma)e dx as. (IT-26) 


The ray is assumed to enter at x = 0. For homogeneous turbulence of a scale A 
with A/5 <1 and with (0n,/0y)? = (0n,/0z)? this yields 


t) 


# Seem ok No(*) \2 : 
a= 40k. 5 {(B5 dx. (111-27) 
0 


(111-26) applies, for example, to the effect of “noise” caused by wind tunnel 
turbulence upon the performance and ultimate resolution of a schlieren system 
and similar problems with electromagnetic waves or sound. 


Optical Measurement of Turbulence 


Kovasznay [26] has shown that a shadowgraph taken of a field of turbulent 
density fluctuations can be used to determine the correlation function of the 
density fluctuations. 

Assume a shadowgraph is taken at x = 0 of a layer extending from « = 0 
to x — I. The relative light intensity distribution h(y, 2) on the plate is related 
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to the refractive index 1 by 
1 


h(y, 2) = const ye pindx, (IIT-28) 
; 


where 


KOovASZNAyY prepares two identical photographic slides. From the transparency 
of these two plates placed face to face he obtains the correlation function h,,- 
by displacing one with respect to the other. From (III-28) then follows, 


I 
hae = h(y, 2) hy +y,2+0) = const / / [V3()] eye W3%e yi nere Ox 46. 
7 0 
(III-29) 


For small scale isotropic and homogeneous turbulence Kovasznay finds, for 
example, that the three-dimensional spectrum of the refractive index fluctua- 
tions E(k) is related to the measured correlation function from the photographs 
by 

k2 E(k) = — | 1 Blr) Jol r) ar ? (III-30) 


0 


where B(7) = hz, — hoo , and Jy is the Bessel function of the first kind and 
of zero order. 


Dipole Scattering of Electromagnetic Waves [21] 


The dipole moment z of a volume v due to an electric field strength F and 
a fluctuation in the dielectric constant D is given by 


x(t) = | D(&, t) Elé, t) de. 


The scattered radiation becomes proportional to 


inl? = | / D(x, t) D(&’, t) E(%, t) E(x’, t) dx dx’. (III-31) 


z* can thus be expressed in terms of the correlation function of D and related 
to turbulent density fluctuations. Specifically, for fluctuations with a length 
scale A small compared with the wave length of the radiation and which varies 
slowly in time we can put 


D(x, t) D(x’, t) = D? A8 d(x — x’) (III-32) 
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and hence 
0 = E? A®y D? (II1-33) 


equivalent to Rayleigh scattering from very small scatterers. 

In more general terms, scattering due to turbulent fluctuations of the 
refractive index can be put into the form of (III-6) by considering the Maxwell 
equations with a refractive index involving a small time and space dependent 
perturbation. The perturbation terms are then considered the forcing function 
(III-6). 

Thus, again writing , for the fluctuation in refractive index we have, for 
example, the equation for E 


= Oz : Ot = 
pep. oe a2 mS ppm). (ls) 


Writing bs E, + @ where ¢ denotes a small perturbation due to the presence 
of the fluctuations 7,, (II1I-34) can be brought in the form of (III-6) 


Lex) = fi (I 11-35) 
with é 
ites ON Be Oe 
i= sang? NM — aire (Eo 7%) 


and the general formalism outlined above applies to the approximate system 
(III-35). 


Production of Sound by Turbulence 


Turbulent fluctuations in an incompressible fluid represent pure shear 
motion or transversal waves. In a compressible fluid the longitudinal pressure 
waves and the transversal shear waves are coupled and hence a turbulent 
field will be coupled with a sound field. The feeding of energy from turbulence 
into sound is a very important problem of turbulence in high-speed and 
especially supersonic flows. Indeed, probably the only important difference of 
turbulent motion at high speeds as compared to low-speed turbulence is this 
coupling with a longitudinal pressure field. The problem of sound production 
from turbulent fluctuations has recently been attacked in a general way by 
LIGHTHILL [30]. In this problem the forcing function in the sound wave equa- 
tions comes from the turbulent, Reynolds stresses. The equation discussed by 


LIGHTHILL is 


where o denotes the density of the fluid and a the velocity of sound. t;; is a 
stress tensor, to which the largest contribution comes from the turbulent 


fluctuations. 


ZAMP III/22 


338 Hans W. LIEPMANN ZAMP 


The radiated acoustic power then involves correlations of the derivatives 


of t,;. It is interesting to note that L1GHTHILL is able to make reasonable predic- 


tions for the dependence of the radiative power on velocity of the flow, Mach 
number, etc. without a detailed knowledge of the turbulent stresses. 


(C) The Effect of Turbulence Upon the Motion of Solid Bodies 


A falling rain drop, an airplane moving in gusty air, smoke stacks in the 
wind, etc. execute on a macroscopic scale a type of motion analogous to 
Brownian motion in small scale. There are, however, essential differences 
mainly due to frequency characteristics of the forcing function as compared to 
the proper frequencies of the systems. In Brownian motion and related pheno- 
mena we mostly assume a white power spectrum of the fluctuating force. A 
white spectrum is equivalent to a 6-function for the correlation function. It 
has been seen above how much the formalism is simplified if the correlation 
length or time is short compared to the characteristic dimensions of the system. 
Another great simplification in Brownian motions results from the fact that the 
probability distribution can be assumed normal or simply related to a normal 
distribution. In turbulence this is often not even a good approximation. Finally, 
the forcing function and the damping in Brownian motion result from the same 
general molecular mechanism. Any dissipation appears in heat and thus in the 
forcing function. This fact is material in the trend toward equipartion of energy. 
In the turbulent counterpart to Brownian motion the forcing function and the 


damping may involve quite different mechanisms. The dissipated kinetic — 


energy of a fluctuating rigid system may go directly into heat or may pass 
through a stage of turbulence. Only the latter case is closely similar to Brownian 
motion. 

Consider, for example, an airplane moving through gusty air as a system 
under the influence of a random forcing function. Then it is clear that the 
system is sensitive mainly to a medium range of frequencies of gusts. The 
very slow gusts are counteracted by the pilot while the very rapid gusts are 
averaged out by the size of the plane. Hence, a problem like this corresponds 
to the behavior of a system, which has a generalized admittance which vanishes 
for very small and very high wave numbers, in a field of forces characterized 
by a specific generalized power spectrum. So far the state of turbulence in the 
atmosphere has not been known with sufficient accuracy and in sufficient 
detail to provide the necessary information for treating such a problem. Much 
progress has been made in the last few years. In some recent investigations the 
response of an airplane was actually used to obtain the spectrum of turbulence 
in the atmosphere [24]. Measurements like these correspond roughly to, say, 
KappLer’s determination of the Avogadro number from the Brownian motion 
of a torsional balance. 


defined by 
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In the example of the airplane it is obviously of great importance to know 
also the probability distribution of the atmospheric turbulence as well as the 
‘filtered’ distribution function of, say, the stress at some part of the structure 
resulting from the gusts. Even if the systems considered are assumed linear 
this problem is by no means trivial. For structural fatigue considerations we 
require also the number of reversals of a given load, a problem related to the 
number of zeros or extrema of a random function. 


(D) The Effect of Turbulence of a Conducting Fluid on a M agnetic Field 


Isotropic, homogeneous turbulence in a conducting fluid and its influence 
on a magnetic field H has been the subject of several recent studies ([19] and 
[23]). In these studies the approximate equations governing the behavior of a 
magnetic field in a moving conductor are usually put into a form similar to the 
KArmAn-HowarTH equation in turbulence (Section V). The statistical proper- 
ties are then discussed along the lines of isotropic turbulence theory. Here a 
different problem is discussed. The approximate equations for the magnetic 
field H in a turbulent conducting fluid with conductivity o and permeability 
are 

0H 


7 A 2H = curl (4 x H) (III-37) 


here A = (4 wo)-! and u(x, t) is the turbulent velocity field. H is assumed 


zero for ¢ < 0. At ¢=0 a constant field H, is switched on. The perturbation . 
field h = H — H, due to the turbulence will be studied. Here we deal with a 
typical example for which the operator includes a time derivative and where 
a homogeneous forcing function is switched on at ¢ = 0. Choosing Hy || x1, we 
have, since div u = 0 


Oh; e. Ou; ts = 
DE Oy ye poet = [LH (I11-38) 


From (III-17) we obtain the following equation for the spectral tensor of h, H,; 


t 
Hy(k, ) = H3 | f dodo! T(k, t— 0) P*(b, to!) Falk, a — 0"). (11-39) 


0 


I'(k, t) is the space Fourier transform of the fundamental solution (7, #) 


SY _ A yty = 82) 6), yO) =0- (III-40) 


Thus 
LR th= wae. (III-41) 
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(III-39) becomes 


Hy(k,t) = Hy | | dodo! e*8-°- Blk, o — 0’) (III-42) 
0 
or 
—Ak*t ; = 
H,,(k, t) = 2 H2 <, i / sinh (4 k? 0) F,,(k, t — 0) do. (I1I-43) 


0 


F’,; is the space Fourier transform of the space-time correlation /,;;. The space- 
time correlation tensor of 0u;/0x, does not appear to be easy to write down in 
general. Thus, only two limiting simple cases are considered here. If 


F,(k, t) = F(R) T dt (111-44) 
then 
H,; = H3 F(R) a [la eneas (ITI-45) 
If 
F,,(k, t) = FY) (III-46) 
then 
H,; = HBF) sep LP. (III-47) 


These two extreme cases differ in the characteristic time; in (III-45) the time 
scale of the forcing function enters while in (III-47) the only characteristic 
time is A k*. Both are proportional to F,,(k) which is related to the spectrum 
tensor @, (hk) of isotropic homogeneous turbulence (see Section IV) by 


—- => 


Fi;(k) = kh} ©,;(R) . (IIT-48) 


The magnetic field increases with time and the perturbation reaches a steady 
state provided of course that the turbulence is kept up by an external agent. 
The ultimate magnetic energy in the perturbation and the spectral distribution 
depends upon the space and time correlation of the turbulence. 

In the cases studied by BATCHELOR [19] and CHANDRASEKHAR [23] the 
increase and steady state of a weak random field is discussed. The comparison 
of their results with (III-42), etc. is therefore not directly possible. The mech- 
anism responsible for the increase in field strength is, however, the same. 
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Impulsive Motion of an Infinite Plate 
in a Viscous Compressible Fluid’) 


By Mitton D. vAN Dyke, Pasadena, Calif.?) 


1. Introduction 


One of the standard solutions of incompressible viscous flow is the motion 
due to an infinite flat plate impulsively started moving in its own plane with 
uniform velocity. Admitting compressibility greatly complicates this simple 
problem. Viscous dissipation heats the fluid, producing pressure and density 
gradients and velocities normal to the plate which culminate in a receding 
shock wave. 

One method of attacking this problem is to linearize the equations of motion 
for small values of the Mach number M. This approach, initiated by LaGER- 
STROM, COLE, and TRILLING [1]’), was generalized and carried to fruition by 
Howartu [2], who succeeded in calculating values of the pressure on the plate. 
The method yields a perturbation field proportional to M?, and perhaps re- 
presents the first step of a series solution in powers of M?. HowartH also 
investigated the initial stages of the motion using the full equations. 

Another method of attacking the problem is to start with the boundary 
layer approximation, which is useful for large time. The resulting solution is 
not limited to small M2, but is valid for arbitrarily high Mach numbers provided 
the time elapsed since the start of the motion is correspondingly large. ILLING- 
WORTH [3] has given the boundary layer solution, and STEWARTSON [4] has also 
briefly discussed its effect upon the outer flow field. 

In the present paper the boundary layer solution is improved in a system- 
atic manner. The procedure involves alternately iterating upon the boundary 
layer solution in the boundary layer and upon the acoustic solution in the 
outer flow field. The two solutions are matched at the outer edge of the bound- 
ary layer by comparing their asymptotic expansions. For simplicity, the as- 
sumptions are made of Prandtl number equal to 1, linear variation of viscosity 


1) Research prepared under Office of Naval Research Contract N 6 ONR-244, Task Order 
VIII, Project NR 061-036. 

2) Guggenheim Aeronautical Laboratory, California Institute of Technology. At present, 
Aeronautical Research Scientist, N.A.C.A. Ames Aeronautical Laboratory, Moffett Field, Cali- 


fornia. 
3) Numbers in brackets refer to References, page 353. 
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with enthalpy, and an insulated plate. The results complement HowarTuH’s 
solution, and agree with it in the region of common validity. The third ap- 
proximation is carried out in detail. Interaction between the boundary layer 
and the outer flow field is found to reduce the skin friction. 

The author is indebted to P. A. LAGERSTROM, who suggested the problem 
and guided the present study, and to Lester LEEs, who offered helpful 
criticism. 


2. Transformation of Equations 


Suppose an infinite insulated flat plate lies initially at rest on the plane 
y = 0 in a viscous compressible heat-conducting fluid. At time ¢ = 0 it is set 
impulsively into motion in its own plane, moving in the positive x-direction 
with uniform velocity U. The velocity components u and v, pressure , density 
g, and enthalpy 7 depend only upon y and ¢. The equations of continuity, 
momentum, energy, and state are 


ort (2 v)y=0, (1a) 

Q (ue + V Uy) = (UW ty)y , (1b) 

Q(%4 +r) => (u%)y— by, (1c) 

elt viy) — (bet op) =n (ue+>03)+2 (Wi)y, (0d) 
p= 2" oi. (Le) 


Here the viscosity « and thermal conductivity 4 depend only upon the enthalpy, 
and the Prandtl number o = ju c,/A is assumed to be constant. The appropriate 
boundary conditions are 


u(y,0)=v(y,0)=0, x4(0,t)=U, (¢ > 0) 
0, 


i(y,0) =i, v(0,2) = 
p(y, 0) = po, i,(0,t) =0, | (2) 
oly, 0) = Qo; 


where subscript 0 indicates values in the initially undisturbed fluid. 
It is convenient to transform the co-ordinate normal to the plate, intro- 
ducing in place of y and ¢ new independent variables Y and T defined by 


| e dy, let. (3) 
0 
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Then 


wt 20 Q 
Vr pe ros Y, , Vy a= en 
Tage as 0 0 0 Qo 0 
Ot or Fut fey men (wee OR ing yf Oye (4) 
Us Maite es Maas 
OG bg n0.ve = Oe Ona O80) 


The role of this transformation is to simplify the operator 0/0¢ + v(0/0y) which 
appears in equations (1). Thus, integrating the continuity equation (la) be- 
tween 0 and y gives 


= et dy = -Y, 
20 : : 20 y : (>) 
so that 
0 Omer 0 
Give oie sou. 9) 


The solution is further simplified by assuming that the viscosity is propor- 
tional to the enthalpy so that according to the equation of state (le) 


ee Py 
Lo Ho Po” “) 


Then the momentum and energy equations (1b) to (1d) reduce to 


tig Ve (2 eile: (8a) 
rain (Zn) 8 0 
i, PF = yb (+H) + 2 (Sotyly: (8c) 


where 1%) = Mo/@o- 


3. Boundary Layer Solution 


For purposes of iterating upon the boundary layer solution, the last three 
equations are written as 


Ei 9, (2 vy), —8r (9a) 


Uz —VoUyy = Vo are : iby) > (9b) 
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45. — Vy UR, eee lyy = % eo uy + cr Vo = V5 
Pr Mf Pe 
+ oi tom ( rm ty), (9c) 


(The first two have been interchanged so that the iteration scheme proceeds 
in order.) 
For large time, that is, large values of the ‘““Reynolds number’’, 
UAt 


VY 


R= 


(10) 


the right-hand sides of these equations can be neglected. Then according to the 
first equation (9a) the pressure is constant through the boundary layer: 


p= bo- ‘ (11) | 


The second equation (9b) is identical with the equation for the incompressible 
problem, which has the solution 


u = U ERFC n (12a) 


where 
y 
og is b 
1 = Vn (12b) 


The general solution of the third equation (9c) for enthalpy has been given 
by ILLInGwortH [3]. However, the result is greatly simplified if the Prandtl 
number o is assumed to be 1. Then, as first shown by Crocco [5], the enthalpy 
is a quadratic function of « alone. With the present boundary conditions 


i=ig+Uu—s a, (13) 


Because the pressure is constant to a first approximation, it follows that 


‘4 = 7 = 7 = 1+ (y—1) M* prec n (1 —— eREC 7). (14) 
Here the Mach number M is the ratio of the speed of the plate U to the speed 
of sound in the undisturbed fluid, which is given by c2 = (y — 1) ig. 

These results are expressed in terms of the distorted normal co-ordinate Y, 
whose relationship to the physical co-ordinate y can now be determined. 
According to equation (4), yy = 09/0, so that integrating equation (14) with 
respect to Y gives 


y=¥[1+(y-]) M? (exec 9 — > ERKc? y 


7A one Lee 
++ jin % eRF 2 = ip" ERE n)|- (15) 
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The »~::ni. velocity v in the boundary layer can now be determined from 
the cortiuuity equation. According to equations (4) and (5), v = yz, so that 
differs niiating equation (15) gives 


ge = eee 2 Vo al =e ; 
ai 1E M ee fz ERF 2 H—€ "ERE 7). (16) 


As usual in boundary layer theory, the condition that v vanishes far from the 
plate has been given up. Infinite values of 7 should be interpreted as referring 
not to the undisturbed flow, but to conditions at the outer edge of the boundary 
layer. Therefore, the normal velocity just outside the boundary layer is 


v | 


eux eet eae s|| ole, 
at 2 ue Vas 0) 


and the displacement thickness 6 is 


)= [thew tt ae 1) M2V» T , (18) 


4. Outer Flow Field 


The velocity component w parallel to the plate vanishes exponentially as 
n > co, so that at the outer edge of the boundary layer there exists only a 
normal velocity component. Therefore, to a first approximation, the flow field 
outside the boundary layer is that which would result from moving the plate 
normal to itself with the velocity given by equation (17). If the structure of 
the receding shock wave were of interest, this problem should be solved in- 
cluding the effects of viscosity. For example, linearizing the equations would 
lead to a piston problem of the sort treated in Ref. [1], Section 2.2. However, 
for purposes of improving the solution in the boundary layer, it is sufficient 
in the first step to treat the outer flow as an acoustical problem. 

For purposes of iterating upon the acoustic solution, equations (1) are 


written as 


Or + 09 Vy = (00 — 2) %y — 4 Oy» 
Oo He = (Qo — Q) Me — QU My + (Mt My) y 


4 
Qo % + py = (Oo 0) % QU + Z(H %y)y> 
: ipa 4. 4 ; 
Oo 41 — Bt = (00 — @) % oi, topy tm (ue + 5 oF) + Hide. 
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Neglecting the right-hand sides leads to the acoustic equation for v: 


—s, : Vet — 0 . (20) 


c 
The appropriate initial and boundary conditions are 


v(y, 0) = v,(y, 0) =0, | 


ioe | ie (21) 
en Pn re Lee 

0(0,2) = TM | NE eer: | 

As usual, the boundary condition has been imposed at the plane y = 0 rather 
than at y = 6(¢). Thus the edge of the boundary layer is represented by 7 = co 
for the boundary layer solution, but by v = 0 for the outer flow solution. The 
solution of this problem is 


eins | mene 
eae 22) 


fort > y/cy = 0 and v = 0 otherwise. The corresponding result for the pressure 
is found to be 


ee , 
P=Pot VEE M? 0 co fget . (23) 


In order that this solution may later be matched with the boundary layer 
solution, it is expanded in series for small values of y/cy ¢, and then expressed 
in terms of the boundary layer variables T and 7, which gives 


y— 1 “Vo 1 j 
Meo W ttm a M? 0904 |/78 Ji+ >. ae 
2 


V Cot 
VT | ly, y—1 
Bat ae ME Gos Yah tg ea pa ld 


5. Second Approximation for Pressure in Boundary Layer 


Having found a first approximation for the outer flow field, we return to 
the boundary layer. There a second approximation can be obtained by iteration. 
Thus, substituting the first approximations for p [equation (11)] and v [equation 
(16)] into the right-hand side of equation (9a) and integrating with respect 
to Y gives 

i Y= 1 aga fore (Sve * BRE n + 


7a : Fi neRK V2 ") +h(T) (25) 
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where A(T) is an arbitrary function which must be determined by matching 
with the outer flow solution. Now for large 7 this is approximately 


Ct Ce # AKL), (26) 


which can be matched with the outer flow solution of equation (24) by proper 
choice of the function (7). The agreement between wipe of order yp is 
deceptive; the result is actually correct only to order »j/*, so that the main 
term in equation (25) cannot be used at this stage. Hence the second approxi- 
mation for pressure in the boundary layer is 


y— 1 VY, 
ee (27) 


At the surface of the plate (7 = 0) the pressure increase is 
Po 2a Tt £9 sa We bs ra 


This is precisely the asymptotic result for large time found by Howartx [2] 
using the method of linearization for small M?. His solution is based upon the 
assumption that o = 3/4 rather than 1, which suggests that the result is in- 
dependent of Prandtl number. 


6. Second Approximation for u in Boundary Layer 


Let the second approximation for tangential recy be «+u’, where w is 
the first approximation given by equation (12), and w’ is a small correction 
term. Substituting the first approximation for w and the second approximation 
for p into the right-hand side of equation (9b) gives 


y if <= 1) — 72 
LT ee Oe Ve ae, M* Vr — URE (29) 


The correction must not disturb the already correct initial and boundary 
values, so that 


u'(y, 0) = u(y, 0) = (0, t) = 0. (30) 

The solution of this problem is found to be 
hes V2 a s|/% Ae ete 3] 
w= -——y(y- I Me ne, (31) 


which again vanishes exponentially as 7 > oc. 
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In the same way, iterating upon equation (9c) leads to a second-order 


correction to enthalpy in the boundary layer to be added to the first approxi- 


mation of equation (13): 
(1 = i M? y e~" ERF n) ; (32) 


Here the initial condition 7 = 7) at f = 0 has been given up in favor of requiring 
the solution to approach the acoustic value at the outer edge of the boundary 
layer. 

Next, improved values of @ and Y can be calculated, and then a second- 
order correction to normal velocity is found which is to be added to the first 
approximation of equation (16): 

- Vi enV mat 


ht eS Bs Fee 
V2 & TF 


? 


xn f + y (vy — 1) M? |pREC y (1 — - ERFC n) at = n e7"* ERF n|| : (33) 


It will be seen later that at the surface of the plate the corrections to u and 
4 compensate each other, so that to second order the skin friction is unchanged. 


7. Second Approximation for Outer Flow Field 


Returning to the outer flow field, we seek a second approximation by 
iteration. Substituting the acoustic solution into the right-hand sides of equa- 
tions (19) leads to a nonhomogeneous wave equation to be solved for v. As 
boundary condition near the plate we require that the solution when expressed 


in terms of boundary layer variables shall for large 7 agree to second order 


(order v9) with the boundary layer solution of equations (16) and (33). This 


matching can be accomplished, which is some assurance of the correctness of _ 


the procedure. This leads to the remarkable result that to second as well as 
first order the outer flow field is just that produced in an inviscid fluid by a 
solid piston moving according to the displacement thickness of equation (18) 
Thus to second order in the outer flow field 


v= 7 =) yrl/ oie. ME satu ae i Lane 
V2a t— y/[Cq ( ), 


2 7 Co t — yey 4 tree yc 
(34) 
va ‘Vo 
= oa Vs = 

p Pot V2a Lo t— y/cy | 
: (35) 

ce ot ): Wie Yo (1 Via Westie 2 yilcy | 

tn 80 Fyfe TOTS e Peon ) 
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8. Third Approximations for p and u in Boundary Layer 


The expression for pressure in the boundary layer given by equation (25) 
remains valid to third order. The arbitrary function f(T) can now be deter- 
mined to this order by expressing the outer flow solution of equation (35) in 
terms of the boundary layer variables and comparing with equation (26). 
Again the matching is successful, with the result that in the boundary layer to 
third order 


P= for TMU ee 


yee! Vo ik eee i} ; (yeh 1) (y= A), 
ik 2 yy a 2 2 
ee Va 00 7 [a (3° Ply eye ERF | n) + arat me). 
(36) 
The pressure increase at the surface of the plate is 
pare er ly Ape? A et ty Se Me 
> ——— a = = Bil 
Po learn F Bx ee ) et 


which is to this approximation just the pressure on the surface of a solid piston 
moving with the displacement thickness. 

Now let uw” be a third-order increment in tangential velocity, to be added to 
the first- and second-order contributions of equations (12) and (31). Substi- 
tuting the second approximation for w and the third approximation for p into 
the right-hand side of equation (9b) gives 


7] tT e (yv oe 1) Me Yo 
Uy: == ANOS i i aa Qa : cs Te 


e™ [Is (1 — 47?) e-" ERF y 


1 < > 2 AOD et = 2) | 
Seg 27) ene V2.5 + 3a " a q(llg es 1 yee 


This is to be solved together with zero initial and Raters conditions. Some 
computation yields the solution 


Pe 38 2 acre rao tewe 2 
tY 22 ree vast V2 (1 2) € ERF J 7 


a 2 1 0) 
+ ne ERE? y= ae" ERE 


4 
7 re 8 2 a 
= (LER V3 = Paar i e° ERF V3 sds) 
a7 | n—2n 


0 


| 


of YM? {By nee” + Va (y +1) (RF 4 —2ne™ [e* eee s ds)}). 
4y/n fs (39) 
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This vanishes exponentially at the outer edge of the boundary layer, but only 
as 73 e~"’, whereas the first and second approximations vanish as 7-1 e~” and 
ye” respectively. Thus each successive term in the boundary layer solution 
penetrates farther into the outer flow field, which is apparently the only way 
in which tangential velocities will arise in this procedure. 


9. Third Approximation for Skin Friction 


The skin friction coefficient is given by 


Sate 25 a) 


The term (f/fp).,~9 18s obtained from equation. (37), and the last term can be 
calculated from equations (12), (31) and (39). The result is 


2 1 rs(y—1 M4 2 
yf sean a 
The first term is the usual boundary layer result. It is independent of com- 
pressibility because of the assumption that viscosity is proportional to tempera- 
ture. Consequently, the decrease of skin friction with Mach number (represented 
by the latter terms) is entirely the result of interaction between the boundary 
layer and outer flow field. 


10. Discussion 


The above procedure can, in principle, be continued indefinitely. In the 
next cycle, a shock discontinuity must be inserted between the outer flow 
field and the undisturbed fluid beyond, and account must be taken of disturban- 
ces reflected inward from the shock wave. Viscous terms in equations (19) will 
become effective at this stage, so that the outer flow field is not inviscid beyond 
the second approximation. At some later stage, the effects of viscosity within 
the shock wave must be considered. 

Clearly the process yields only an asymptotic solution for large time. It 
is not valid near t= 0, because equation (41) would imply that an infinite 
impulse per unit area is required to set the plate in motion. 

The boundary layer approximation is permissible only if the maximum 


normal velocity is small compared with U. According to equation (17), this 
requires that 


7 a1 (42a) 
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Likewise, acoustic theory is applicable to the outer flow only if the normal 


velocity at the edge of the boundary layer is small compared with the speed of 
sound, which requires that 


= Zen. (42b) 


Hence the present solution is a perturbation scheme for small M 2/VR or 


M3] V R, whichever is the larger. The preceding results suggest that the solution 
can be expanded in powers of these two parameters. 
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Zusammenfassung 


Eine unendlich ausgedehnte Platte werde zur Zeit ¢= 0 in ihrer eigenen 
Ebene mit konstanter Geschwindigkeit in Bewegung gesetzt. Das betreffende 
Strémungsfeld einer zahen, kompressiblen Fliissigkeit wird untersucht. Die 
Lésung geht von der Grenzschichtnaherung aus und benutzt ein Verfahren sukzes- 
siver Annaherungen, um den Effekt der Druckwelle, die wegen der Kompressi- 
pilitat entsteht, darzustellen. In zweiter Naherung ergibt sich eine verbesserte 
Grenzschichtlésung und eine akustische Lésung ausserhalb der Grenzschicht. 
Die in der erstgenannten Lésung auftretende unbestimmte Funktion Jasst sich 
aus der Forderung berechnen, dass beide Loésungen am Rande der Grenzschicht 
iibereinstimmen sollen. Aus der Wechselwirkung der dusseren Stromung mit der 
Grenzschicht ergibt sich eine Verminderung der Schubspannung an der Platte. 
Das benutzte Verfahren ist nur fiir grosse Werte eines dimensionslosen Zeit- 
parameters giiltig. 
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La plaque rectangulaire fléchie d’épaisseur linéairement variable 


Par HENRY Favre?) et BERNHARD GILG, Zurich 


Différents auteurs, notamment R. GRAN OLssoN, ont étudié la flexion des 
plaques vectangulaires d’épaisseur variable*). Le cas ot cette variation est 
linéaire n’avait cependant pas encore été traité jusqu’ici. Or une telle variation 
peut étre facilement réalisée, surtout dans le béton armé, en choisissant pour 
les faces de la plaque deux plans formant un petit angle. Nous avons donc estimé 
utile d’étudier ce cas. 

Aprés avoir établi l’équation différentielle régissant la flexion (§ 1), nous 
indiquerons deux méthodes pour intégrer cette équation lorsque la plaque est 
appuyée le long du contour (§2). Nous ferons ensuite une application de ces 
méthodes, en supposant que la surcharge soit une pression hydrostatique (§ 3). 
Dans un dernier paragraphe (4), nous donnerons sous forme de diagrammes 
les résultats d’un calcul numérique. 


§ 1. Equation différentielle d’une plaque rectangulaire 
d'épaisseur linéairement variable 


Soient a et b les cétés. Choisissons les axes x, y, z indiqués dans la figure 1 
et posons, pour l’épaisseur # en un point quelconque: 


=|1+4 (72 -1)| ho, (1) 
d’ou | 
(yo =(L—A)hy, (B)y.p= (+A) h 


Ay désignant la valeur de l’épaisseur de long de la droite y = b/2. Le nombre A 
est la variation, divisée par hy, que subit h depuis cette droite jusqu’a l'un ou 


1) Chaire de Mécanique, Ecole polytechnique fédérale. 

*) Citons entre autres: R. GRAN OLsson, Die Biegung der rechteckigen Platte verdnderlicher Dicke, 
Thése de l’Ecole polytechnique de Berlin (1933); Biegung der Rechteckplatte bei linear verdnderlicher 
Biegungssteijigkeit, Ing.-Arch. 5, 363 (19384); Biegung der Rechteckplatte von exponentiell verander- 
licher Dicke, Bauingenieur 21, 230 (1940); Biegung der Rechteckplatte bei linear verdnderlicher 
Steifigkeit und beliebiger Belastung, Bauingenieur 22, 10 (1941), — E. Reissner, Remark on the 
Theory of Bending of Plates of Variable Thickness, J. Math. Phys. 14, 43 (1937). 
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Vautre des cétés qui lui sont paralléles. Nous supposons donc que / varie 
uniquement dans la direction de l’axe y. 


Z. 
> 
of a ae the 
pet 
0) £, ae GO 
feutllet moyen 
| 
: 
yy 
Fig. 1 


L’équation régissant la flexion d’une plaque ot l’épaisseur varie selon une 
loi quelconque est, en coordonnées cartésiennes rectangulaires!) : 


: oD oO ODI Owl eee 
DAAC + 25+ ep AS +25, ° Gy AS + AD AS | 
, TDD LAS ee lel Bal age > cae caa | 2) 
ae) as "Oy? Ox dy Ox dy | Oy? sa) =P, 
ou 
ee E h?® 3) 


1Z(1 297) ? 
£ étant le déplacement d’un point du feuillet moyen, A, l’opérateur de LAPLACE 
02/dx2 + 0?/0y?, E, le module d’élasticité, v, le nombre de Potsson, et P(x, ¥), 
la surcharge par unité de surface. 


1) Voir S. TimosHENKO, Theory of Plates and Shells (McGraw-Hill Book Company, New-York 
et Londres, 1940), p. 195. 


356 Henry Favre et BERNHARD GILG ZAMP 


Dans le cas étudié, l’équation (2) devient, compte tenu des relations (1), 
(3) et en posant pour abréger 


_ 332 E 14¢ (1 — ») a + ; (4) 
+ (ee = 1)" AAt] a (=; ats 1) Fa (4c — (1— >») 5 
i a (42 1) + At + (= 23 1)’ AAG| : | 


Cette équation différentielle n’est évidemment applicable que sil’épaisseur / 
est partout petite par rapport aux autres dimensions et si le déplacement ¢ 
d’un point quelconque du feuillet moyen est lui-méme petit par rapport a h, 
ce que nous supposerons essentiellement dans la suite. Nous admettons égale- 
ment qu’aucune force n’agisse dans le plan de ce feuillet. 

Les expressions des moments de flexion M,, M, et de torsion M,,, par 
unité de longueur, sont les mémes que lorsque |’épaisseur est constante: 


z hel RES uel ee Hes eee Me 
M,=—D (s+ ¥553)) My = Dee re at) | 
Myy=—(—»D as on 
Il en est de méme de celles des tensions oy, 0, T,, en fonction des moments: 


M, ° M 
On = 2, % = 


Sree eri v ats M gy 
etme ka 


Wyig*?. tay = haya, 7 (6) 
z désignant la distance d’un point de la plaque au feuillet moyen. En particulier, 


les valeurs de ces tensions aux points z = h/2 de la face opposée 7 a celle ot 
agit la surcharge # sont données par les formules: 


6M. 6M 6M 
(2); = h2 =, (%)i= he ", (Try)e = ae” ae (7) 


§2. Méthodes d’intégration, dans le cas ou la plaque est appuyée le 
long du contour 


Le déplacement ¢ dépend non seulement de x, y, mais aussi de A. Dévelop- 
pons ¢ en une série illimitée, procédant selon les puissances croissantes de cette 
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_derniére variable: 


C(x, y, A) = Colx, vy) + Ci(%, vy) A+ G(x, y) A? +e = oS Cah. (8) 


1=0 


Pour déterminer les fonctions ¢,(x, y), introduisons l’expression (8) dans les 
deux membres de l’équation (4), puis identifions les coefficients des mémes 
puissances de 4. Nous obtenons le systéme d’équations différentielles simul- 
tanées: 


we| AAt, = £, | 
fa| AAC, = —3 [5-5 Abo + (4% — 1) 446), 
3] Adt, = —3[F- Ge 40 + (FX - 1) 446] 

~3 [Fe {to- 1 -») Geet + (4% — 1) gy At 


East ciieris ta Mehie: titel felts; el aiteqces a eOiel else io; ol felie. o) cel omesief (eo ey(e) 16° <a" 'e) tes jeh ie) to) 8) 0) ©, Bye eye eye 8 9, 'e 


(9) 
#| AAC; = 3 |5 is Ati. + (27-1) Adee] 
~3[ 2 (aces—a—») Sesh + 5 (25 - 1) 9g Abs 


La plaque étant appuyée le long du contour, les conditions aux limites sont: 
C=0 ae = 3 =0, (10) 


n désignant la direction de la normale au contour, au point considéré. Pour 
que ces conditions soient satisfaites quel que soit A, il est necessaire, d’aprés (8), 
que I’on ait en tout point du contour: 


G0 ete = y= Ly Ze see) (we) 
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Revenons au systéme (9). La premiére équation contient seulement la 
fonction inconnue C9. Cette équation n’étant autre que celle d’une plaque 
rectangulaire posée d’épaisseur constante ji), nous pouvons Vintégrer par la 
méthode de Navier. Dans ce but, commencons par développer le second membre, 
qui est une fonction connue de x, y, en une série double de sinus, en posant: 


pee “, a ee Mae . nay 
ia Hx, v) Pays, Amn Sin——— sin eet (12) 
ou 
a b 
4 . MeE* ‘ = 592. 76.9) 
AS wae PW is sin —~ ax / f(x, v) sin a dy. (13) 
0 0 


Posons ensuite: 


CO [oe) 
at) EFC es nmUuY 
Ce pts ay oe SL errata Some ane (14) 


m=1 n=1 


En introduisant les séries (12) et (14) dans la premiére des équations (9), et 
en identifiant, nous obtenons pour les coefficients: 


|: i ia Amn 2\2° (15) 


a‘ D, ex | 7; i: 

On vérifie facilement que la valeur ainsi obtenue pour ¢, satisfait aux con- 
ditions aux limites (11). 

En introduisant cette valeur dans le membre de droite de la seconde des 
équations (9), ce membre devient une fonction connue de x, y. On peut alors 
le développer en une série de sinus et déterminer ¢,, comme on l’a fait pour Cy. 
En introduisant ensuite, dans le membre de droite de la troisiéme des équations 
(9), les valeurs obtenues pour ¢, et ¢,, on est ramené a des opérations analogues 
aux précédentes, pour le calcul de fy, et ainsi de suite. 

La méthode est tout a fait générale. Son application suppose seulement que 
la fonction p(x, y) soit développable en série, ce qui est toujours le cas dans la 
pratique. Cependant, une difficulté se présente dans les calculs. Si Cy S'exprime 
a l’aide d’une série double, les fonctions Cy, cy, ... se présenteront en réalité 
sous la forme de sommes de séries, dont une seule est double, les autres ayant 
un ordre de multiplicité supérieur & deux. L’exemple traité dans le §3 le 
montrera éloquemment. 

Les séries multiples étant relativement peu commodes pour les calculs 
numériques, il sera donc indiqué, toutes les fois qu'on pourra le faire, de chercher 
a appliquer la méthode de Maurice Livy, Celle-ci consistera ici 4 développer les 
seconds membres en séries simples et 4 mettre également les fonctions ¢; sous 
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la méme forme. Cette seconde méthode ne sera toutefois applicable que si la 
surcharge p a des valeurs égales en deux points symétriques par rapport a la 
Broite. «= 2/2. 

L’exemple traité au paragraphe suivant, ot il est possible d’appliquer les 
deux méthodes, fera bien ressortir leur différence. 


§3. Cas d’une pression hydrostatique 


Supposons que deux cétés soient verticaux et que la plaque limite un 
bassin entiérement rempli d’un liquide de poids spécifique y (figure 2). On a, 
dans ce cas: 


a he (16) 
(Mh 
— . —- sf 


K-La 
(1+A) ho 
| 


Fig. 2 yV 


Premiére méthode 


La premiére des équations (9) devient: 


AAE=cy, =“ 
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en posant, pour simplifier, 
Peet (18) 


En remplagant, dans la formule (13), f(x, y) par c y, et en intégrant, on 
obtient: 


8 b c m = Lhe Ee sh * ? co 
_ fe] ya 
a ( 1) eta , t= ie eH ae ee) 
d’ou, en vertu de (14), (15) et en posant 
g= +: (19) 
es 8 arc = = (— 1)"-2 - Mae . ny 
o= a n sin > (2a 


La seconde des équations (9) s’écrit, compte tenu de (17): 


Wat 3 eee ye teat 21 

emia Far ote (4% — ) 9]. (21) 
En partant de l’expression (20) de Co, on obtient sans difficulté: 

0 _ 8arc va = (— 1)” MO Mae 3 nary 


Développons cos (7 x y/b) en série simple de sinus, en appliquant des for- 
mules analogues 4 (12) et (13): 


nny — nay 
COS — = . A., sin eee: 


b : " be 
n=) 
ou 
bo aI cos { x (n’ + n)} 1 — cos{ 2 (n’— n)} ; 
A | | a gt ale eee 


< | _ 1 cos{x (n’+ n)} 
a (n’ + n) ; 


Introduisons ce développement dans (22), ce qui donne pour OAE,/Oy la 
série triple : 


fi) 8 2 fo, 3) [o,2] lo, 2) ney n 
Te tae ee PD 


Sis 8 a 
m=1,8,.. w= 1,2, .. 9’ Ire, % (m* + n? 9?) 


n+ n n’—n a b 
eee 


1 — cos{ 2 (n’ + n) 1— cos{a(n’— n)}] . max . ns 
x fe OR: j igh Meni s NS Lakh Sin 2S tei SEN 
=0,sin’=n 
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En développant (2 y/b — 1) y en série double, on obtient: 
2 _ 166 7 ol | 
SP PaO eee | 


(—1)” 2 en UOMEEE 0 TOE 
x [= "= = cos (na) Seer cosa) | Si sin 


(24) 


Le membre de droite de |’équation (21) devient ainsi, en tenant compte de 
(23) et (24), la somme d’une série triple et d’une série double. Posons en consé- 
quence: 


co foe) 0O : 
? E DY ; Da mann 210 ne, oan sin — ome 
= . ne V 
2 De Das See 


Substituons les expressions (23), (24) et (25) dans l’équation (21), puis iden- 
tifions. Nous obtenons directement les valeurs des coefficients C,,,4 ¢t Dinn 
qui, introduites dans (25), donnent finalement, aprés quelques transformations: 


BATE CM 802 nes oS = (— 1)” 
es SDD ame 


m=1,3,... n=1,2,.. n’=1,2,... 


(25) 


1 il so EET GED 
x ; — ; sin sin 
ni+n n—n a b (26) 
=0, p. val. paires des dénom. 
Ail if 1) 
5 = ae nn? ne 5D TREE. EE 
> A sin Sit— A 
+= pe mn (m® + n® 9?)? a b | 


Cette fonction satisfait aux conditions aux limites (11). 

Nous avons également calculé ¢,, qui est la somme d’une série quadruple, 
d’une série triple et d’une série double. Il n’y aurait aucune difficulté a déter- 
miner ¢, et méme ¢,. Ces calculs seraient toutefois moins utiles, car, dans la 
pratique, l’exactitude obtenue en limitant le développement (8) aux trois 
premiers termes suffira en général. 


Deuxieme méthode 


Posons: 
Coe meron SUN el (27) 
W=13 5008 
ou 
a, = Cees (28) 
a 
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Yon désignant une fonction de la variable y. La série (27) satisfait aux conditions 
(11) le long des cétés x = 0 et x = a. Introduisons cette série dans la premiére 
des équations (9), en multipliant en outre le second membre par 


oo 


1= = ee = sin (a, x); mm V45)) 
nous obtenons: 
oo oc 
Dy (WE — 208 Yo, + of Y,) sin (a, 4) = 2% 5" = sin (a, 2), 
n=1,38,... n=1,3,... 


ou Yo, et Yjy désignent la dérivée seconde et la dérivée quatriéme de Y,,, par 
Tapport a y. 
En identifiant, dans les deux membres, les coefficients de sin(%, *), on est 
conduit a l’équation différentielle ordinaire: 
J > 4cy 
Yon — 20%; You + at Yon = Ga 
La solution générale de cette équation est: 


Ven a ‘Aon cosh (Xp y) 5 Bon sinh (&n y) i Con Ln y cosh (mn y) | 


30) 
A 4cy ( 
Dy n&y V sinh (a, y) 46 a a3 ? | 
ou Ay», ..., Do, désignent des constantes. 
Les conditions aux limites (11) exigent que 
ton 0 St ky, = (31) 


pour y=Oet y= 0b). 
Ces conditions, exprimées a l’aide de la solution (30), conduisent 4 quatre 


a x? sinh («,, b) ’ 


équations linéaires qui, résolues par rapport aux constantes Aon, ..., Don 
donnent les valeurs: 
4bec ba 
Ay, = 0 Io. == . a 
an » Bon a «> sinh («,, b) (1 3 tgh («,, d) : | (32) 
2be¢ 
Con = Don = 0. | 


En introduisant ces valeurs dans le second membre de (30), puis en substi- 
tuant dans (27), on obtient finalement : 


ee Ste ome pagegos EA (+53. 
A ao |b sinh («,, 6) 2 tgh (a, d) 


= 1 Bis 


(33) 
CAR VECOSE: (io 0) ine | 
2 sinh (a, 6) | sin (On x) . 
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Posons ensuite: 


j= Y VERS cea ea i (34) 
7 — os 
ou 
nN I 
OO rae (28) 


Y,, ne dépendant que de y. La série (34) satisfait aux conditions (11) le long 
mes-euies=7 — 0 et 4 ="a. 

En introduisant les séries (27) et (34) dans la seconde des équations (9), 
celle-ci devient, compte tenu des relations (17) et (29): 


D, Win — 205 Mi + o%,,).sin @,%) 


n=1,3,... 
4a! 4 4c¢ (29? 
=— z a, 5 E Con cosh (a y) — = al aaa ( an — sin (%, %) 


ou C,, est la constante donnée par la troisiéme des formules (32). 
En identifiant, on obtient pour Y,,, l’équation différentielle ordinaire: 


VEG, aly, 


ln n 


a -3 (7% [2 E ieuaiu togny = ae |+ oe (rz = yj}, 


6 
a a, 


dont la solution générale est: 


Yin = Ayn cosh (% ¥) + By» sinh (% ¥) + Cin %n y cosh (on ¥) | 


8s) 
+ Dyn %n y sinh (4p, y) + G+ Gy y* cosh (oy y) + Gay + Ga y*, 
oul 
— 48¢ 6¢ DelZ2 eee 
Tg ae re ora a a4 sinh (a, 6) ’ Gs ee Cae abob’ (36) 
Ain +++ Dyn tant des constantes. 
Les conditions aux limites exigent que 
Y= Oy .cureva, = 0 (37) 
pour y = 0et y=). Onen déduit : 
48 ¢ Sod shia SSRIS 9 Bee) 
Ain = aba,’ Baer a «, sinh (a, b) | b a2 fs , 38) 
6b¢ ee oor | 
Cin = a a, sinh (a, 6) ’ Die a& sinh (x, b) * 
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1 et conséquem- 
Les fonctions Y,, sont donc poeta ere déterminées, q 


ment ¢). 
Posons encore 
~ 
Ce oo Yon sin (Xn x) , (39) 
n=1,3,... 
ou 
On =, (28) 


Y,, ne dépendant que de y. ' 
En procédant comme nous l’avons fait pour Y), et Y,,, on obtient pour les 
fonctions Y,,: 


Xan = Aan cosh (tn.¥) ati sBen SIND (ety) Te Came Y Ai a | 
+ Dyn %& y sinh (a, y) + H,+ Hy + A, y?+ A, v8 (40) 
+ H, y® cosh (%, y) + H, y? sinh (%» y) + Hz y* cosh (a ¥) , | 
ou 
Hm ree he | | 
A= Gea e~ apr: Ee pee aaeee ne 41) 
f= a ERY (- Fa : i mate Ay ahah cay 
A435 = aber Pen = aes ah by? 
Bin =~ sarmecy [37-39 + SE a + ea 
96 cosh (ay b) 3(7— 119) . 18 (eee 6 1, Liga 
BF a2 ~ ba, tgh (a, b)  tgh® (a, b) tgh® (a, b) 
SO ear Saha “a tH v3 a8 OE) : 
Seaaead | 


C2 est ainsi connu. Pour la raison invoquée a propos de la premiére méthode, 
nous n’avons pas calculé d’autres fonctions ¢;. 

En introduisant les valeurs (20), (26), ... ou (33), (34), (39), ... de Cg, C, «+ 
dans la formule (8), on obtient le déplacement ¢ en fonction de x, y. Les for- 
mules (5) permettent ensuite de calculer les moments de flexion et de torsion 
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M,, M, et M,,, Vou l’on déduit les tensions a l’aide des relations (6) et (7). 
Nous ne donnerons pas ici les valeurs explicites de ces grandeurs, pour ne pas 
alourdir notre exposé. On pourrait d’ailleurs calculer également les réactions 
réparties le long des cétés, ainsi que les réactions concentrées agissant aux 
quatre sommets du rectangle. 

Remarque. Nous avons supposé que la surface libre du liquide était au méme 
niveau que le c6té supérieur de la plaque (figure 2). Les deux méthodes que 
nous avons indiquées permettent cependant de résoudre sans difficulté le cas 
ou la surface libre est située au-dessous de ce cdté, et celui ot elle est située au- 
dessus. 

Ce dernier cas est évidemment le plus simple, car il suffit d’ajouter, a la 
valeur de ¢ que nous venons de calculer, celle correspondant a une pression 
constante fy, agissant sur toute la surface de la plaque. On trouve par exemple, 
en utilisant la seconde méthode et en posant 


Powe 
Da = Co. (43) 
ae, | on cosh ( (Xn y) oP Bon sinh (On y) ai Con An V cosh (%n y) | 
. 4 : 
+ Donn y sinh (%n ¥) + 5 “| Sil ares | 
ou 
4¢ -. == I ‘ 2 Co 
Aon mains ot ar , Bon a ‘sinh (a), b) af, a [1 cosh (On b)| | és 
fy b C 26g. l= cosh (%p, 6) e 2e, | 
4 ee 2 sinh ( mk on Gia sinh (%,0)- ’<"°" — a, a” 
Cy = AN [Ain cosh (%n y) a ss sinh ee y) ar Cy n&n V cosh (On y) | 
w=1,3,.. ; (46) 
+ Dyn %m y sinh (%, ¥) + G, + Gy y + Gy? cosh (an y) 
+ G, y? sinh («, ¥)] sin (%» x) , 
ou 
IAG 6¢, 1+ 2cosh (a, 5) 
Ain = — a, ri poi ot, ae sinh («,, b) ‘ 
6 Co ; hia. b sinh a2) | 
Cin = G5 a sinh («, 6) 1 — cosh (an 8) + “a7 an 
6¢ Sea _ 1265 Cpe 24 Cy 
Hse Ee E a, b sinh (%» 0) o ae ae ai ab’ 
6)  1— cosh (Xp 8) ee AO Gor 
Gs ogo | sinha, b). ~ Ce at ab 
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Co= D7 [Asn cosh (mY) + Bon sinh (on ¥) + Con &m ¥ COSh (on 9) 
N= 1,3, 


+ Don %m y sinh (a, y) + Hy + Hy y + Hz v?+ Hy y? cosh (a, Y) 


+ Hy y? sinh (a, y) + Hg vy? cosh («, y) + H, y3 sinh (a y)| sin (a, x) , 


ou 
24¢ 4 (1+ ») 
Asn = — aie 1+ area 
Pais Co { 4»— 8+ 2 (7 — 5 ») cosh (a,,-d) 
Ben x5 a sinh (a, b) | %” | sinh ta, 0) 
12 (1 — ») [1 — cosh (a, b)] cosh? b 
af sinh? (a, b) es | ples agate 
I COshi(o 0) a2 36 (7 + 3 ») [1 — cosh (a, b)] 
pote eS | sinh (a, 6) | zs “a, 6 sinh (a, 6) 
a 96 (1 — ») [1 — cosh (a, b)]) 
25 ie 
an) Oe | emery oe 
Co Daa chta wee 36 — 12 vy + 24 cosh (a, 8) 


(1 — ») [1 — cosh (a, b)] 1 
sinh? (a, b) a | 


ae a“ |4 sinh (x, byes 4 (23 + 5 ») [1 — cosh (Xp =] 


sinh (a, 6) 
+ art (7 + 3 v) [1 — cosh (a, yy}. 
De Es fi mene). 305538) 
ma Boe h5M), a5, 
fig aot pe Nd ae -at |" ace ee a Ra : - " : 
Os ee wat 4 es Tae 3 


bi 10 (5 + ») [1 — cosh (a, d)] 
%, b sinh («,, b) | . 
2 €) (7 + ¥) [1 — cosh (a, d)] 2 7 
Hae 0 n _ 26)(7+ ») 
. ay, @ 6? sinh (a, b) » = ms ab 
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§ 4. Application numérique 


Revenons au cas de la figure 2. Nous avons fait les calculs numé€riques 
complets, par la seconde méthode, en supposant la plaque carrée (a = b) et 
A = 0,21). Pour le nombre de Poisson, nous avons pris y = 0,25. 

Les diagrammes de la figure 3 représentent les valeurs obtenues pour les 
grandeurs ¢, M, et M,, le long de l’axe vertical de la plaque (x = a/2), pour 
trois approximations du calcul. Dans la premiére, on a limité la série (8) a 
son premier terme, dans la deuxiéme, A son second terme, etc. 

On voit que la convergence est assez rapide et qu’il serait inutile de calculer 
une quatriéme approximation. C’est pourquoi tous les calculs ont été faits en 
définitive en ne gardant que les trois premiers termes de la série (8). 

La figure 4 montre les courbes de niveau des surfaces représentant les fonc- 
tions ¢, M,, M, et M,,. Il est intéressant de constater que la surface ¢ est quasi 
symeétrique par rapport au plan horizontal y = 0/2, ce qui montre que la plaque 
étudiée se déforme sensiblement comme si son épaisseur et sa surcharge étaient 
constantes. Cela signifie que la variation d’épaisseur compense, en quelque 
sorte, la variation de la pression, quant aux déformations. 

Les surfaces M,, M, et M,,, sont, par contre, loin de présenter ce caractére. 
Ce fait est d’ailleurs compréhensible, car la partie inférieure de la plaque, étant 
moins déformable que la partie supérieure, accuse des moments plus grands 
que dans cette derniére, ott les pressions hydrostatiques sont d’ailleurs plus 
faibles. : 

Enfin, la figure 5 montre les courbes de niveau des surfaces (0,),;, (dy); et 
(Tz); - Les deux premiéres surfaces sont, comme celle de ¢ (figure 4), sensible- 
ment symétriques par rapport au plan horizontal y = b/2. Cela est di a un 
effet de compensation, analogue 4 celui que nous avons signalé plus haut a 
propos de ¢. La surface (t,,); n’accuse malheureusement pas ce caractére de 
symétrie: dans la partie inférieure de la plaque régnent des tensions tangen- 
tielles beaucoup plus grandes que dans la partie supérieure”). 

Les caractéres de symétrie que nous avons pu relever, qui existent malgré 
l’asymétrie de la surcharge par rapport a l’axe y = 6/2, montrent éloquemment 
le role joué par la variation linéaire d’épaisseur de la plaque. Ce réle est favo- 
rable, puisqu’il contribue a donner une bonne répartition aux déformations et 
aux tensions normales dans le domaine considéré. 


1) Ils’agit, par exemple, d’une plaque dont l’épaisseur varie de 8 a 12 cm. 
2) Les tensions tangentielles ne sont pas nulles le long du contour. Ce fait entraine, comme 
dans le cas ot l’épaisseur est constante, l’existence de réactions supplémentaires qu'il serait facile 


de calculer. 
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£7, 
Surface a 


Si 
a x 
f -> 


Surface. <—£ 
f ry x 


Fig. 4 
Plaque de la figure 2. Cas ot. a = b, A = 0,2 et vy = 0,25. Courbes de niveau des surfaces ¢, M,, M, 
a g 2. 3 
et M,.,,. 
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Fig. 5 


Plaque de la figure 2. Cac aod 5 
Plaque de la figure 2. Cas ot a = b, A = 0,2 et » = 0,25. 


Courbes de niveau des surfaces : 
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Summary 


This paper deals with bending of simply supported rectangular plates of 
linearly varying thickness. The solution of the partial differential equation of 
the middle surface may be reduced to the integration of a system of simpler 
equations by development in a series. Any of these equations has a form similar 
to the equation of a rectangular plate of constant thickness. The solution may 
be found either by the method of multiple trigonometrical series (NAVIER) or, 
for symmetrically distributed load, by simple series (MAURICE LEvy). 

The case of hydrostatic pressure is discussed in detail. A numerical example, 
with graphical results, clearly shows the influence of variability of thickness on 
bending, on twisting and bending moments, and on normal and shearing forces. 


(Regu le 13 juin 1952.) 


Die automatische Kursregelung von Flugzeugen 


Von Rupotr Graus, Minneapolis, Minn.*) 


(Fortsetzung und Schluss) 


7. Die automatische Kursregelung eines Flugzeuges in der Langsachse 


Bei einem Langstreckenflug wird ein Flugzeug meistens in einer konstanten 
Hohe geflogen. Dem Piloten stehen Instrumente, wie der Héhenmesser und 
das Variometer, zur Verfiigung, welche ihm Anderungen in der vertikalen 
Flugrichtung anzeigen. Wenn in konstanter Hohe geflogen wird, dann ist der 
Anstellwinkel « gleich dem Winkel %. Hierbei ist # definiert als der Winkel 
zwischen Langs- und Horizontalachse des Flugzeuges, wahrend der Anstell- 
winkel « als der Winkel zwischen Anstrémrichtung und Flugzeuglangsachse 
bezeichnet wird, wie aus Figur 7 ersichtlich. 

Infolge Brennstoffverbrauches andert sich der Anstellwinkel « bei konstant 
gehaltener Motorenleistung. Damit das Flugzeug auf gleicher Hohe bleibt, muss 
@ gedndert werden. Aufgabe der automatischen Kursregelung ist es jetzt, einen 
kommandierten Winkel #, so einzuhalten, dass sich die Hohe des Flugzeuges 
nicht andert. Hierbei kann der Eingang 0, vom Piloten selber oder aber auto- 
matisch durch einen Héhenregler verandert werden. Der Ausgangswinkel 0, 
wird durch einen Vertikalkreisel gemessen, wie in Figur 7 angedeutet. Das 
Blockdiagramm der ganzen Langssteuerung ist in Figur 8 dargestellt. Die Ab- 
weichung «, ist wiederum die Differenz zwischen #, und #,. Sie wirkt tiber 
einen Regler auf eine Verstirker-Servomotor-Kombination. Der Ausgang des 
Servomotors ist iiber Kabel mit dem Héhensteuer verbunden. Ein Héhen- 
steuerausschlag 6, ergibt durch das Flugzeug eine Anderung von #,, die, wie 


+) Minneapolis-Honeywell Regulator Company, Aeronautical ‘Division. 


Ww 
~I 
bo 


Rupo.r GLAaus ZAMP 


A 


Fig. 7 
Schematische Darstellung eines Flugzeuges in der Langsachse. 


VERTIKAL 
KREISEL 


Fig. 8 


Blockdiagramm fiir Flugzeuglangssteuerung. 


bereits erwahnt, vom Vertikalkreisel gemessen wird. Damit erhalten wir einen 
geschlossenen Regelkreis. Alle Komponenten im Blockdiagramm von Figur 8 
sollen als gegeben betrachtet werden, mit Ausnahme des Verstarkungsfaktors K p 
des Reglers. Dieser soll so eingestellt werden, dass das maximale Amplitudenver- 
haltnis | 0/9.) (¢ ) | gleich 1,3 wird. Es ist ebenfalls nach der Eigenfrequenz w, 
des geschlossenen Regelkreises gefragt. Als erstes werden jetzt die Gangkurven 
der einzelnen Regelelemente bestimmt werden. 


Regler, Verstdrker und Servomotor 


In unserm Beispiel sind sowohl Regler wie Verstarker und Servo elektrische 
Komponenten des C-1-Autopiloten, hergestellt durch die Minneapolis-Honey- 
well Regulator Co, Die Reglergangkurve ist vollstandig beschrieben durch einen 
variablen Verstarkungsfaktor Kp, da die dynamischen Eigenschaften vernach- 
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lassigt werden kénnen. Der Verstarkerfrequenzgang ist ebenfalls nur durch eine 
Konstante bestimmt. Hingegen weist der Servo dynamische Eigenschaften auf, 
die nicht vernachlassigt werden diirfen. Die Ausgangsseite des Servo treibt iiber 
ein Reduktionsgetriebe die Seiltrommel und ein Potentiometer. Am Potentio- 
meter wird eine Spannung proportional der Seiltrommeldrehung abgegriffen 
und zum Verstarkereingang zuriickgefiihrt. Somit wird jeder Eingangsspan- 
nung des Verstarkers ein bestimmter Drehwinkel der Seiltrommel zugeordnet. 
Der Frequenzgang 6;/e; der Verstarker-Servomotor-Kombination wurde ex- 
perimentell bestimmt und ist in Figur 9 dargestellt. 


T Fa tiraet ob 
GANGKURVE gs (iw) DES_SERVOMOTORS 


° Ss 
Es GE VERSTARKUNGSFAKTOR Kg = Odb IN GRAD/ GRAD 


AMPLITUDE |6s| IN Dezibel 
E 


9,1 10 KREISFREQUENZ w IN SEC” 


Fig. 9 


Das Blockdiagramm in Figur 8 kann jetzt vereinfacht werden, indem 0s/és 
als ein Block, wie die gestrichelte Linie andeutet, gezeichnet wird. 


Das Hohensteuer 


Das Héhensteuer wird mittels Drahtseilen durch den Servomotor bewegt. 
Als Flugzeug soll in unserem Beispiel eine DC-3 verwendet werden. Da die 
Eigenfrequenz von Hohensteuer und Drahtseilen viel hoher ist als der Frequenz- 
bereich, in welchem die iibrigen Komponenten arbeiten, wird das dynamische 
Verhalten dieser Teile vernachlassigt. Hingegen ist zu beachten, dass der H6- 
hensteuerausschlag fiir einen konstanten Drehwinkel der Servoseiltrommel mit 
der Fluggeschwindigkeit variiert. Der Steuerausschlag Ist am gréssten am 
Boden und nimmt mit zunehmender Fluggeschwindigkeit infolge der auf das 
Héhensteuer einwirkenden Luftkrafte und der Elastizitat der Drahtseile ab. 
Da wir in unserem Fall nur einen Flugzustand betrachten, reduziert sich die 
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Gangkurve des Héhensteuers zu einem konstanten Verstarkungsfaktor 


K, = 0,44 Hohensteuerausschlag in Grad 2 Sp Dodatvele (27) 


‘Servodrehwinkel in Grad 


Das Flugzeug 


Die dynamischen Eigenschaften eines Flugzeuges sind abhangig von Flug- | 
geschwindigkeit, Hohe, Gewicht und Schwerpunktslage. Deshalb muss eine 
optimale Regeleinstellung entsprechend diesen variablen Parametern geandert 
werden. Wir betrachten aber hier nur den Fall, wo Fluggeschwindigkeit, Hohe, 
Gewicht und Schwerpunktslage als konstant vorausgesetzt sind. Somit kann 
das Verhalten des Flugzeuges in seiner Langsachse durch eine Gangkurve be- 
schrieben werden. Der Frequenzgang Ky Gp wurde im Flugzeug direkt experi- 
mentell durch sinusférmiges Variieren des Héhensteuers aufgenommen und ist 
in Figur 10 aufgezeichnet. Er kann ebenfalls aus der Differentialgleichung des 


dC -3 - 
GANGK URVEN Siw) DES FLUGZEUGES 


FLUGGESCHWINDIGKEIT v= 230 km” - 


HOHE = 2200 m UBER MEER 
GEWICHT = 10100 kg 


VERSTARRUNGSFAKTOR Ke = [3 db IN GRAD/GRAO 
AMPLITUDE |. 
ri et 


Ss) 
| PHASENVERSHIEBUNG «p 
a ES IN GRAD 


T0,0 


KREISFREQUEN2 W IN SEC” 


Fig. 10 


Flugzeuges abgeleitet werden, falls die aerodynamischen Konstanten fiir diesen 
bestimmten Flugzustand bekannt sind [9], [10]"). Der analytische Ausdruck fiir 
die Flugzeuggangkurve Ky F rip) hat die Form: 


ps 0 ao + ap-+ a, p? 
Ky Gry = —* () = — a. 
FF we 28 
(e) Oe () by + b, p + by p? + by p? +- by pt (28) 
1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 382, 


Vol. III, 1952 Die automatische Kursregelung von Flugzeugen Sys) 


Hierbei sind die Werte fiir die verschiedenen a und 6 durch die aerodynamischen 
Konstanten bestimmt. Die Berechnung des Frequenzganges mit Gleichung (28) 
ist sehr miihsam, da die Ausdriicke fiir die darin auftretenden Konstanten 
einen komplizierten Aufbau haben. Andererseits ist eine direkte Bestimmung 
des Frequenzganges im Flugzeug zeitraubend und kostspielig. Es wird deshalb 
oft die Antwortkurve des Flugzeuges aufgenommen und daraus der Frequenz- 
gang bestimmt. 


Der Vertikalkreisel 


Dieser Kreisel errichtet eine vertikale Bezugsachse im Raum und ist als ein 
Pendel mit einer sehr langen Schwingungsdauer zu betrachten. Die Gangkurve 
des Vertikalkreisels kann durch eine Konstante Ky, beschrieben werden, da 
das uns interessierende Frequenzspektrum viel kiirzere Schwingungsperioden 
hat als der Kreisel. In unserm Beispiel ist: 


Grad 


if — rege — “ a @) \ 
Ko tna = 0 Dezibel. (29) 


Die Gangkurve K Gi) des aufgeschnittenen Regelkreises 


IS ISt. 
FONG en py Foe 


(t@) S(i@) Ky Ky Cri w) Kyx (30) 


und 


ee et Kora Ka eK ys DEI bels (31) 


-. of 
‘Saal 


Hierbei ist Kp noch unbestimmt, wahrend die tibrigen Verstarkungsfaktoren 
gegeben sind. Es wird: 


K=K,+6 in Dezibel, (32) 
Gis) ae Csi) Grin) (33) 
oder . 
| Gee | = | sua ae | Gria) in Dezibel , (34) 
P(Giw)] P[4sGay] + PGR o)]" 3) 


Gio) ist im Nichols-Diagramm, Figur 11, aufgetragen. 

Die Kurve Gi) wird jetzt so lange in der Vertikalrichtung verschoben, 
bis sie Tangente an die M-Kurve 1,3 = 2 db ist. Der Ordinatenunterschied der 
beiden Kurven ist gleich dem Verstarkungsfaktor K des aufgeschnittenen 
Regelkreises und wird aus Figur 11 abgelesen: 


Be a= Om enivel: (36) 
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Fig. 11 o[+88-] IN GRAD 


Gangkurven fiir Flugzeuglangssteuerung. 


daraus kénnen wir den grésstméglichen Verstarkungsfaktor Kp des Reglers 
mit Gleichung (32) bestimmen: 


: Grad 
Kp = —2,5 db = 0,75 Guat (37) 
Die Eigenfrequenz , des geschlossenen Regelkreises wird ebenfalls aus F igur 11 
abgelesen : 


®,= 3,75 5-4, (38) 


In Figur 12 ist die Gangkurve des geschlossenen Regelkreises fiir das Amphi- 
tudenverhaltnis und die Phasenverschiebung tiber den Logarithmus der Fre- 
quenz aufgetragen. Diese Darstellungsart ergibt ein iibersichtliches Bild tiber 
die Regelgiite des geschlossenen Systems. 

Es ist ersichtlich, dass die Regelung iiber einen Frequenzbereich von 0,1 bis 
2,5 s~tein Amplitudenverhaltnis von ungefahr 0 Dezibel und vernachlassigbarer 
Phasenverschiebung hat. Uber diesen Frequenzbereich vermag der Ausgang 
dem Eingang sozusagen verzégerungsfrei zu folgen. Es soll jetzt untersucht 
werden, wie das System in Figur 8 verbessert werden kann. Oder anders aus- 
gedriickt : 

Der Verstarkungsfaktor K und die Eigenfrequenz @, sollen erhéht werden, 
und zwar wiederum fiir ein 
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+ 
GANGKURVEN DER GESCHLOSSENEN REGELKREISE z (iw ) 
FUR SYSTEM I UND : 
SYSTEM I: OHNE WENDEZEIGER —— 


SYSTEM Ir: MIT WENDEZEIGER -—e 


pel 


Se 1,0 KREISFREQUENZ ® IN SEC 10,0 


Die Vergrésserung des Verstaérkungsfaktors K in einem Regelkreis hat eine 
dreifache Wirkung: Erstens wird der Effekt von Stérungen auf das statische 
Verhalten der Regelung verkleinert, zweitens die Geschwindigkeit, mit welcher 
das System bei einer plétzlichen Eingangsanderung antwortet, erhéht, und 
drittens die Stabilitat des Regelkreises vermindert. In unserm Fall hat K 
bereits den maximalen Wert erreicht, der mit unserm Stabilitatskriterium 
toleriert werden kann. Infolgedessen brauchen wir ein zusatzliches Regelele- 
ment, um die Gangkurve des aufgeschnittenen Regelkreises K Gjj) So zu 
formen, dass der Verstarkungsfaktor K erhéht werden kann. Als zusatzliches 
Regelelement verwenden wir einen Wendezeiger, der nichts anderes als ein 
federbelasteter Kreisel ist und die zeitliche Anderung ‘von #, misst. Dadurch 
ergibt sich ein Blockdiagramm, wie aus Figur 13 ersichtlich ist. 


HOHEN LUGZEUG 
STEUER eee 


VERTIKAL 
KREISEL 


Fig. 13 


Blockdiagramm fiir Flugzeuglangssteuerung mit Wendezeiger. 
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Es ist: ; 
joer 8 


Der grundsitzliche Aufbau eines Wendezeigers ist aus Figur 14 ersichtlich. 


w 
n 
x 
oO 
a 
4 
o 
= 
1) 
xc 


FEDER Fig. 14 


Schematische Darstellung eines Wendezeigers. 


Der ideale Wendezeiger verméchte allen zeitlichen Anderungen von @, fiir jede 
beliebige Frequenz zu folgen. Es wiirde dann der Ausgang: 


O, = Ki Pere (40), 


Hierbei ist Ky die sogenannte Kreiselkonstante und in unserm Beispiel als 1 
angenommen. Bei der praktischen Ausfiihrung eines Wendezeigers kénnen die 
dynamischen Eigenschaften von Kreiselmasse, Federn und Reibungskraften in 
der Lagerung des Kardanringes im allgemeinen nicht mehr vernachlassigt wer- 
den. Sie werden mit Hilfe eines Systems zweiter Ordnung in die Gangkurve 


des Wendezeigers einbezogen. Wir erhalten dann: : 
ge = We, 0 U 
i. (Sa) ok Ae Wee Fie (41) 


Es bedeuten: 


w = Stérfrequenz s~! 
Me», — Eigenfrequenz von Wendezeiger s~} 


5 = Verhiltnis von tatsachlicher Dampfung zu 
kritischer Dampfung des Wendezeigers. 
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In unserm Beispiel ist 


Gy, = LOS ancl =a. 


em 


Die Gangkurve (Q,/9,) (( w) wird erhalten durch Ersetzen von w durch Bem 
und ist in Figur 15 aufgetragen. 


i lapels le alee enrol a 
= 
; —|GANGKURVE ° DES WENDEZE IGERS FUR Kz 
a Ke bciu LCN We, = 10 SEC 
ae | (cea = gE .02 4 
eS 
30 Wem +100 


PHASENVERSCHIEBUNG @p 
IN GRAD 


\ 
4+20 
a 
| | a fae 
0,1 1,0 KREISFREQUENZ GIN SEC? 100 
Fig. 15 
Wir erhalten jetzt: 
&y = a; = Oe ) ; (42) 
& =a—Q,, (43) 
f= &s — 05 - (44) 


Alle Gangkurven und deren Verstarkungsfaktoren sind gleich wie im vorigen 
Beispiel, nur dass jetzt die Gangkurve von (Q,/9,) (¢ w) zusatzlich dazukommt. 
Die Bestimmung der Verstarkungsfaktoren Kp und K, der beiden Regelkreise 
erfolgt so, dass die Regelung bei B zuerst aufgeschnitten und der innere Regel- 
kreis, durch Bestimmen von Ky , eingestellt wird. Es wird systematisch so vor- 
gegangen : 
a) Zeichne die aufgeschnittene Gangkurve (Q,/es) (¢ @) ins Nichols-Diagramm 
und stelle den Verstarkungsfaktor Ky so ein, dass M = 2 Dezibel ist. Aus 
Figur 16 erhalten wir: 


F Grad 
Ky = —16,5 db = 0,15 Grad s-}° 49) 


380 Rupotr GLAus ZAMP 


Q KG IN GRAD 
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Gangkurven fiir Flugzeuglangssteuerung mit Wendezeiger. 


b) Berechne (#,/a) (¢ w). Aus dem Nichols-Diagramm kann sofort Q,/a abge- 
lesen werden. Daraus erhalten wir: 


ies ae (eye 
+. (tw) = a (tw) 5 a ae Gs A. (46) 


_ Figur 13 kann jetzt wie in Figur 17 in ein Blockdiagramm vereinfacht werden. 
c) Die Gangkurve des Vertikalkreisels ist gleich 1, und die Gangkurve des 
Reglers besteht ebenfalls nur aus einem Verstarkungsfaktor Kp, der vorerst 
gleich 1 angenommen wird. Damit erhalten wir durch Eintragen von A ins 


REGLER 
Ka 


VERTIKAL~ 
KREISEL 


Kvk 


Fig.'17 


Vereinfachtes Blockdiagramm der Flugzeuglaingssteuerung mit Wendezeiger, 
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Nichols-Diagramm (Figur 16) und Verschieben dieser Kurve, bis sie M = 2 db 
bertihrt, direkt Ap als Ordinatenunterschied. Da A Amplituden aufweist, 
die — 20 db iiberschreiten, ist in Figur 16 als Amplitude (A + 10 db) auf- 
getragen. 


Wi hal 
ir erhalten Grad 


und 


Open io- ok (48) 


Tabelle 2 ist eine Zusammenstellung der Verstarkungsfaktoren K und Eigen- 
frequenzen w, der beiden Regelsysteme. 


smaibeller2 
| Verstarkungsfaktor des auf- Eigenfrequenz des 
geschnittenen Regelkreises geschlossenen Regel- 
aaa =1 
Ran | ie kreises W, S 
System I (ohne Wendezeiger) 25 | iL) 3S 
System II (mit Wendezeiger) 24 | 17 7,0 


Aus Tabelle 2 ist ersichtlich, dass der Verstarkungsfaktor K in System II 
ungefahr zehnmal grésser ist als in System I, wahrend die Eigenfrequenz fast 
_verdoppelt wurde. 

Die Gangkurve des geschlossenen Regelkreises (9,/0,) (1 w) mit Wendezeiger 
ist ebenfalls in Figur 12 eingezeichnet. Aus diesem Diagramm ist ersichtlich, 
dass die Amplitude bei w = 1 s~! von 0 db auf — 4 db abzusinken beginnt, um 
dann auf das Maximum von + 2 db anzusteigen. Dieser Amplitudenverlaut 
wiirde nicht den gewiinschten Giitegrad der Regelung ergeben. Der Grund 
dieser negativen Amplitudenausbuchtung ist darin zu finden, dass die Eigen- 
frequenz des Wendezeigers zu tief gewahlt wurde. Aus Figur 16 ist ersichtlich, 
dass die Gangkurve bei w = 5s~! zu steigen beginnt. In diesem Frequenz- 
bereich macht sich die Resonanz des Wendezeigers bemerkbar. Durch Wahlen 
einer héheren Messwerkeigenfrequenz, zum Beispie! w,,, = 20s~1, kann das 
Ansteigen der Gangkurve herausgezégert werden. Damit kann G;;,,) noch weiter 
nach oben verschoben werden, und die Annaherung an die Kurve M = 0 db 
wird wesentlich verbessert. Mit dieser Anderung wiirde der gewiinschte Giite- 
grad der Regelung erreicht. 


8. Schlussfolgerungen 


a) Die Berechnung einer Regelung und deren optimale Einstellung wird we- 
sentlich erleichtert durch Verwenden der Gangkurve des aufgeschnittenen 


Regelkreises. 
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b) Mit Hilfe des Nichols-Diagramms kann unverziiglich die Gangkurve des ge- 
schlossenen Regelkreises sowie dessen Stabilitat, Eigenfrequenz w, und 


‘ 


4 


maximaler Verstarkungsfaktor K bei gegebenem M-Kriterium bestimmt — 


werden. 

c) Falls eine Regelung den gestellten Anforderungen nicht geniigt, so muss 
die Gangkurve K G durch Zufiigen weiterer Regelelemente umgeformt wer- 
den. Diese Umformung wird wesentlich erleichtert durch das Nichols-Dia- 
gramm, da nur darauf geachtet werden muss, dass die Gangkurve K G iiber 
einen moglichst grossen Frequenzbereich der Kurve M = 0 db folgt. 

Die Praxis hat erwiesen, dass in vielen Fallen mit Hilfe des M-Kriteriums 
direkt von der Gangkurve K G auf das Verhalten der Antwortkurve ge- 
schlossen werden kann. Hiermit wird die Entwerfung und Anderung von 
Regelsystemen, bei vorgeschriebener Regelgiite, wesentlich vereinfacht, da 
nur die Gangkurve des aufgeschnittenen Regelkreises bendtigt wird. 


a1 
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Summary 


This paper describes a method for the analysis and optimum design of servo- 
mechanisms. The transferfunctions of the open loop as well as the Nichols’ dia- 
gram are used for the systematic and rapid solution of control problems. As an 
example of the application of this method, the calculation of an automatic control 
for the longitudinal axis of an aircraft is carried out. The simplest case is first 
considered, using a vertical gyro to control the pitchaxis of the airplane. In 
order to improve this system, a rategyro is added. This more complex system 
serves at the same time as a sample for the investigation of a multiloop system. 
(Eingegangen: 19. November 1951.) 
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Kritertum der Gleitbereitschaft metallischer Vielkristalle 


Von F. ROHNER, Neuhausen am Rheinfall?) 


1. Einleitung 


Wenn bei der Gleitung Bindungen zwischen den betroffenen Atomen zu 
lésen waren, miissten ideale Metallkristalle Elastizitatsgrenzen in der Gréssen- 
ordnung von 1000 kg/mm? aufweisen. Die Elektronentheorie der Metalle hat 
nun aber gezeigt, dass die metallische Bindung als kollektive Wechselwirkung 
zwischen dem negativen Elektronengas und den positiven Atomriimpfen zu- 
stande kommt, dass also gar keine lokalisierten Bindungen zwischen einzelnen 
bestimmten Metallatomen existieren. Arbeit gegen die Bindungskrafte miisste 
bei der Gleitung deshalb nur geleistet werden, wenn diese zwangsweise mit 
einer bleibenden oder auch nur voriibergehenden Volumzunahme des gleitenden 
Metallkristalls, etwa einer Dilatation in der Normalen zur Gleitebene, verkniipft 
ware. Das ist aber nicht der Fall; je nach der Lage des Gleitsystems zur Bean- 
spruchungsrichtung kann die Gleitung sogar mit einer voriibergehenden Normal- 
kontraktion verlaufen?). Die Erklarung des Widerstandes der Metallkristalle 
gegen Gleitung bzw. die Ableitung der Grdésse ihrer Elastizitatsgrenze muss 
deshalb auf anderer Basis gesucht werden. 

W. L. Brace’) hat 1942 die Ansicht ausgesprochen und begriindet, dass 
ein Kristallit erst dann gleiten werde, wenn er durch den Gleitprozess seine 
elastische Spannungsenergie verringern kénne. Die elastischen Spannungen 
werden auf den Kristalliten durch seine Nachbarkristallite tibertragen; dies 
gilt sowohl fiir innere wie fiir 4ussere Spannungen. BRAGG fasst diese einem 
Einzelkristalliten durch seine Nachbarn aufgedriickten Spannungen als Schub- 
spannungen auf und leitet ab, dass Gleitung eintreten k6énne, sobald in einem 
Kristalliten die Schubspannung o den der Gleichung (1) entsprechenden Betrag 
erreiche: 


Gas. (1) 


Dabei bedeuten G, s und L den Gleitmodul bzw. den Atomabstand bzw. den 
Durchmesser des Gleitelements in der Gleitrichtung. 

Rouner?’) fiihrte 1947 aus, dass — da sich der Kristallit bei Zugbeanspru- 
chung im elastischen Bereich als untergeordneter Bestandteil seines Kollektivs, 


1) Forschungsinstitut der Aluminium-Industrie-AG. 

2) E. Scumip und W. Boas, Kristallplastizitat (Springer, Berlin 1935), S. 284. 
3) W. L. Brace, Nature 149, 511 (1942). 

4) F, Rouner, J. Inst. Mct. 73, 303 (1947). 
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des Vielkristalls, betrage — der technische Elastizitaétsmodul des Vielkristalls 
und nicht der Gleitmodul des Kristalliten in die Formel eingehen miisse. 
Schliesslich sei noch auf die Arbeit von Woop und RACHINGER verwiesen}) ; 
diese Autoren schliessen sich den Anschauungen von BracG an. Bei allen 
bisherigen Ableitungen wurden Querkontraktionen, die sowohl bei der elasti- 
schen Dehnung wie bei der plastischen Gleitung auftreten, vernachlissigt. 
Zweck der vorliegenden Arbeit ist es, das Kriterium der Gleitung oder, besser 
ausgedriickt, das Kriterium der Gleitbereitschaft unter Einbeziehung dieser 
Querkontraktionen zu prazisieren. 


2. Ableitung der Ansdtze zur Berechnung der Spannungsenergien 


Um unnotige Komplikationen zu vermeiden, nehmen wir an, dass wir es 
mit einem aus einigermassen isometrischen Kristalliten aufgebauten Vielkristall 
zu tun haben und dass — wie etwa bei Aluminium — keine betrachtliche Aniso- 
tropie der elastischen Eigenschaften in den Einzelkristalliten besteht. Zur wei- 
tern Vereinfachung fiihren wir unsere Uberlegungen an einem wiirfelférmigen 
Kristalliten durch, dessen Querschnitt durch Figur 1 dargestellt sei. Wir denken 
uns diesen Kristalliten als Bestandteil seines Kollektivs, des Vielkristalls, einer 
eindimensionalen Dehnung in dem durch einen Pfeil angedeuteten Sinne unter- 
zogen. Bei einer bestimmten kritischen Dehnung wird eine Gleitung, etwa 
entlang der eingezeichneten Gleitebene, erfolgen. Hernach wiirde der Kristallit, 
falls er nicht zwischen seine Nachbarkristallite eingespannt ware, die in Figur 2 
dargestellte Form aufweisen. Allgemein ausgedriickt: Wenn die Dimensionen 
des wiirfelf6rmigen Kristalliten der Figur 1 


C= Lay 1b Sera 


sind, so sind die Dimensionen des Kristalliten nach der Gleitung gemiiss Figur 2 


Ss 
yous 


e=1+ 7, b=1 a=1, 
wobei s und L Atomabstand bzw. Durchmesser des Kristalliten oder Gleit-} 
elements in der Gleitrichtung bedeuten. Wenn die Gleitung bei der kritischen 
Dehnung « erfolgt ist, muss der Zustand Figur 2, damit die Passung mit den 
Nachbarkristalliten (die auf Grund ungiinstigerer Orientierung und Grosse erst 
spater gleiten sollen) erhalten bleibt, durch elastische Gegenspannungen auf 
folgende Dimensionen zuriickgefiihrt werden: 
Oe Le Bee eee nee ee 


,’ 


m m~ 


Das sind die Dimensionen, die der Kristallit vor Einsetzen der Gleitung — im 
elastisch gedehnten Zustand — aufwies. Dabei bedeutet m die Poissonsche Zahl. 


') W. A. Woop und W. A. RAcHINGER, J. Inst. Met. 75, 571 (1949). 
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Fig. 1 Fig. 2 
Querschnitt durch wiirfelf6rmiges Gleit- Gleitelement von Figur 1 nach der Gleitung. 
element mit eingezeichneter Gleitebene. 


Die Einpassungsverformung wird sich zwar nur zum Teil an dem verformten 
Kristalliten, zum andern Teil aber an seinen Nachbarkristalliten auswirken. Uns 
interessieren hier jedoch nur die energetischen Verhaltnisse. Fiir deren Ableitung 
k6nnen wir uns denken, die Einpassungsverformung erfolge voll am verformten 
Kristalliten. Effektiv werden Gleitung und Einpassung gleichzeitig miteinander 
erfolgen. Gleitung wird nur moglich sein, wenn die zur Einpassung aufzuwen- 
dende Verformungsenergie nicht grésser ist als die zur Erzielung der Dehnung ¢ 
in den Kristalliten hineingesteckte Energie. Unsere Aufgabe besteht also in der 
Ermittlung der Verformungsenergie fiir die Zuriickfiihrung des spannungslos 
gedachten Kristalliten nach der Gleitung (Tabelle 1, Zustand III) auf die 
Dimensionen des Kristalliten im elastisch gedehnten Zustand (Tabelle 1, Zu- 
stand II). Diese Verformungsenergie der Einpassung ist dann mit der Verfor- ° 
mungsenergie der Dehnung von Zustand I auf Zustand II zu vergleichen 
(Tabelle 1). 
Tabelle 1 


Dimensionen des Einheitskristalliten im Ausgangszustand I, im elastisch gedehnten Zustand II 
und, nach der Gleitung um einen Atomabstand, im elastisch vollstindig entspannten Zustand III 


Kristallitendimensionen 


Zustand 
c | b a 
I 1 1 1 
II 1l+e 1 — e/m 1 —e/m | 
III 1+ s/L 1—s/L 1 


ZAMP II1/25 
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Die Verformungsenergien berechnen sich nach der Formel (2): 


Eu, =G |e + + + (€, + & + &3)?|. (2) 


m— 2 


E\., bedeutet dabei die Verformungsenergie pro Volumeinheit, G den Gleit- 
modul, der mit dem Elastizitatsmodul F in folgender Beziehung steht: 


Fiir den Spezialfall 


Spey, Co. ES onay Cao == 


(Verformungsenergie beim Ubergang von Zustand I auf Zustand IT) reduziert . 
sich Formel (2) auf Formel (3): 


m+1 


Eq = Get. (3) 


3. Berechnung und Vergleich der Spannungsenergien 


; \ 
Die Dehnungen fiir den Ubergang von Zustand III auf Zustand II (Tabelle 1) 
betragen: 


3 
oy at i ar ee Lye 
Par Ss € Ss (1 = 
a EN m }’ 
€ Sere = 
Co — = 
3 m Lom 


wb} 


Ss 


wobei e* = € 


Wir kénnen bei allen Dehnungen deren gemeinsamen Faktor s/L weglassen: — 


das bedeutet, dass wir als Einheit der Dehnung s/L wahlen. Ferner lassen wir 
aus den Berechnungen der Spannungsenergien den Faktor G weg, was darauf 
hinauslauft, dass wir die Verformungsenergie statt fiir das Volum 1 fiir das 
Volum 1/G erhalten. 

Die Poissonschen Zahlen der meisten Metalle liegen zwischen 3 und 4. 
Tabelle 2 gibt die Lésungen fiir m — 3, Tabelle 3 fiir m = 4. In Figur 3 sind die 
Ergebnisse dieser Berechnungen zusammenfassend dargestellt. Die Spannungs- 
energien fiir die primére Dehnung sind durch strichpunktierte Linien darge- 
stellt, die Spannungsenergien der Einpassung durch voll ausgezogene Linien. 
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Zu jeder Dehnung ¢ ergibt sich auf der strichpunktierten Linie die Dehnungs- 
energie E(. — jm, ¢/m) und auf der durchgehend ausgezogenen Linie die Einpas- 
sungsenergie E(,. .,..,) , die aufgewendet werden miisste, wenn nach Erreichung 
der Dehnung e* Gleitung erfolgen wiirde. Solange die durchgehend ausgezogene 
Linie der Einpassungsenergie iiber der strichpunktierten Linie der Dehnungs- 
energie verlauft, ist Gleitung nicht mdglich. Es ergibt sich somit, dass fiir m = 3 


15 


O5 L0XYL 


Fig. 3 
Verlauf der zur Einpassung eines Kristalliten in seine Umgebung erforderlichen Spannungsenergie E 
als Funktion der Dehnung ¢, bei der Gleitung erfolgt, und im Vergleich zur Verformungsenergie, 
die zur Dehnung ¢€ gehort; solange die Einpassungsenergie grésser ist als die zu € gehorende Verfor- 
mungsenergie, kann Gleitung nicht erfolgen. 
—.—- zur Dehnung ¢ gehérende Verformungsenergie, 
Einpassungsenergie. 


Tabelle 2 
Dehnungsenergien Eve, pers und Einpassungsenergien Ee a ES a m = 3 und ver- 
schiedene kritische Dehnungen e*, bei denen Gleitung supponiert wird 


e* Eve, —e/m, —e/m) et ee ey Eve, 5, €3) 
0 0 —1 it 0 2,00 
0,3 0,12 =0.7" 0,9 —O7% 132 
0,6 0,48 = 04. (0809 = 0,2 0,88 
0,9 1,08 —0,1 0,7 — 0,3 0,68 
EZ, 1,92 OFZ Os — 0,4 0,72 
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Tabelle 3 
Dehnungsenergien Ev, rim Zahm und Einpassungsenergien Bae, be fiir m = 4 und ver- 
schiedene kvitische Dehnungen e*, bei denen Gleitung supponiert wird 


e* Eve, —elm, —e/m) ef es 

0 0) —1 1 

0,3 0,125 —0,7 0,925 —0,075 1,376 
0,4 0,20 —0,6 0,9 — 0,10 1,20 
0,6 0,45 —0,4 0,85 —0,15 0,95 
0,8 0,80 —0,2 0,8 — 0,2 0,80 
1,0 1725 0,0 0,75 —0O,25 0,75 
12 1,80 +0,2 0,70 —0,30 0,80 


- 


= 


; 


Gleitung nicht méglich ist, bevor e den Wert 0,75 s/L erreicht hat. Fiir m = 4 


liegt der entsprechende Wert bei e¢= 0,8 s/L. Damit haben wir einen untern 
Grenzwert fiir die Elastizitatsgrenze gewonnen. Zur Erzielung der Dehnung 
e = 0,75 s/L ist die Spannung o = 0,75 E s/L erforderlich. Da Gleitung bei 
kleinern Spannungen nicht eintreten kann, muss die Elastizitatsgrenze o; 


grosser sein als 0,75 E s/£. Nach unserer Erfahrung, die wir durch Beispiele — 


soeben noch belegen werden, liegt bei kalt bearbeitetem oder an sich feinkérni- 
gem Aluminium die 69 92-Grenze in der Gegend von 2 E s/L, wenn fiir L rént- 
genographisch oder elektronenmikroskopisch bestimmte Kristallitendurchmes- 
ser eingesetzt werden. Bei diesen Messungen bleiben aber immer vorhandene 
intrakristallin eingebaute Fremdatome und Gitterbaufehler, die die freie Weg- 
lange der Gleitung ebenso wie Kristallitengrenzen verkiirzen, unberiicksichtigt. 
Die Lage der 69 :-Grenze des Aluminiums beim zirka 2,5fachen des derart 
berechneten untern Grenzwertes fiir o, erscheint also durchaus plausibel. 

Die Grosse L, definiert als mittlere freie Weglange der Gleitung oder als 
mittlerer Durchmesser der Gleitelemente, ist — wie soeben angedeutet wurde — 
schwer erfassbar und tragt zu einem guten Teil die Schuld daran, dass wenig 
experimentelles Material vorliegt, das zur Nachpriifung dieser Gesetzmissig- 
keiten geeignet ware. Aussere Dimensionen der Proben, Korngréssen, Mosaik- 
blockgréssen, Konzentrationen von Fremdelementen in fester Lésung, Aus- 
scheidungen und Gitterbaufehler haben alle Einfluss auf die Grésse L. An sehr 
feinen Drahten kann ein Uberwiegen des ersten Faktors, der aussern Proben- 
dimensionen, beobachtet werden, was die oft auffallend hohen spezifischen 
Festigkeiten solcher Feindrahte erklart. Bei Reinmetallen sind die letztauf- 
gefiihrten Einfliisse ausgeschaltet oder reduziert. 

Bei Reinaluminium kann durch sehr ausgiebige Kaltbearbeitung die Kri- 
stallitengrésse bis zu 1 Kristallitendurchmesser reduziert werden, Weiter 


— 
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herunter kommt man nicht; wahrscheinlich vergréssern sich kleinere Kristalli- 
ten spontan durch Rekristallisation. Die Pulvermetallurgie erlaubt aber, Rein- 
aluminiumk6rper mit noch kleinern Kristalliten gewissermassen synthetisch 
herzustellent). Deren mittlerer Kristallitendurchmesser kann aus elektronen- 
mikroskopischen Untersuchungen des zur Sinterung verwendeten Pulvers 
ziemlich genau abgeleitet werden. Wir haben fiir zwei solcher pulvermetallur- 
gischer Erzeugnisse aus Reinaluminium, tiber die ndhere Angaben in der er- 
wahnten Arbeit!) zu finden sind, mit 0,5 bzw. 0,3 « mittlerer Kristallitengrésse 
die 69 92-Grenze gemessen. In Tabelle 4 sind die erhaltenen Werte, im Vergleich 
zu maximal kaltbearbeitetem Reinaluminium, den aus der Formel 6992 ~2E s/L 
berechneten gegentibergestellt. Die gemessenen Werte liegen durchgehend etwas 
hoher als die berechneten. Ein Vergleich der drei gemessenen 09 o)-Grenzen 
unter sich zeigt, dass sie sehr gut dem Gesetz gehorchen, nach dem sich die 
Elastizitatsgrenze umgekehrt proportional zu L verandern muss. 


Tabelle 4 
Gemessene und nach dey Formeél 0) 9, ~2£E s/L berechnete Elastizitdtsgrenzen von Rein- 
aluminium verschiedeney Kristallitengrosse L 


(E = 7000 kg/mm?, s = 2,85 A) 


e kg/mm? 
berechnet | gemiessen 
Aluminium 99,5 . . Ih fei 4 5 
SiMe WGI B Roars 5 0,5 8 ili 
SEPA A) 3 eo oe 0,3 wu its} 15 
1) Geschiitzte Markenbezeichnung fiir gesintertes Aluminiumpulver 


Die in dieser Arbeit entwickelte Betrachtungsweise fiihrt also zur Auffas- 
sung, dass der Widerstand der Metalle gegen plastische Verformung nur in- 
direkt, iiber den Elastizitatsmodul, mit den Bindungskraften des Gitters zu- 
sammenhdangt, dass er aber in erster Linie durch die Einpassungsspannungen 
bestimmt wird, die in jedem Gleitelement — wenn der Zusammenhalt mit den 
Nachbarelementen gewahrt bleiben soll — aufgebracht werden miissen. Diese 
Auffassung fiihrt zur Konsequenz, dass ein Idealkristall von makroskopischen 
Dimensionen einer Verformung praktisch keinen Widerstand entgegengesetzt, 
also schon bei minimsten Spannungen zu gleiten beginnt. Diese Konsequenz 
vertragt sich aber gut mit der Erfahrung. Idealkristalle von makroskopischen 
Dimensionen sind allerdings praktisch nicht realisierbar. Man ist aber heute 
mit der Raffination vieler Metalle und der Ziichtung fehlstellenarmer Kristalle 


1) R. SH ee A, vy, ZEERLEDER und F, Rouner, Festschrift RoS (Vogt-Schild AG., Solothurn 
1950), S. 77, 
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so weit, dass einer Extrapolation auf die Eigenschaften des wahren Ideal- 
kristalls in vielen Fallen keine Bedenken mehr entgegenstehen. Man kommt 
dabei immer auf eine verschwindend kleine Festigkeit des Idealkristalls. 


Summary 


The first section of the paper is devoted to a brief discussion of some principles 
of the electron theory of metals, which show that there are no definite localized 
bonds between any definite atoms of the metallic lattice. Slip is possible without 
any bonds being disrupted or even any energy spent against bonding forces. In 
a crystallite making part of a polycrystal slip is possible as soon as the elastic 
strain energy can, by slipping, be reduced. On this assumption and taking into 
consideration elastic and plastic transverse contractions a stress value is derived, 
which can be considered as a lower limit for the advent of slip. It is inversely 
proportional to the diameter of the crystallite or more exactely to the free path 
of slip and to some extent dependent on Poisson’s ratio. Determinations of the 
elastic limit of wrought and of sintered pure aluminium samples are given and 
shown to be in accordance with this theory. 


(Eingegangen: 10. Mai 1952.) 
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Bemerkungen iiber den Einfluss der Plasmaschwingung 
in Laufzeitréhren 


Von Fritz Lupr, Baden?) 


1k 


Es wird gezeigt, dass die in einer fritheren Arbeit iiber die Travelling-Wave- 
Tube*) beniitzten Ausgangsgleichungen folgende Erweiterungen ohne zusdtzliche 
Ableitungen iiber Raumladungskrafte enthalten: 

a) Plasmaschwingung (Plasmafrequenz), 
b) freier Einelektronenstrahl, 

c) Doppelstrahlverstarker, 

d) Travelling-Wave-Tube. 

Dabei kommt als wesentliches Resultat heraus, dass die Plasmaschwingung 
(das heisst die Raumladungskrafte) nur bei b und c, aber nicht bei d einen 
Einfluss hat. Dieses Resultat steht offenbar im Widerspruch mit dem anderer 
Autoren und sei damit zur Diskussion gestellt. 

Die Ausgangsgleichungen der Arbeit [I] seien nochmals ohne das Ersatz- 
schema angeschrieben. 


1) AG. Brown, Boveri & Cie., Baden. 
2) W. Frey und F. Ltp1, ZAMP 1; 237 (1950) (hier zitiert; Arbeit {1]). 
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Fiir die Helix 


CMU eae (1) 
edge x= 50), (2) 
avi dye 0. (3) 
Fiir den Elektronenstrahl 
OVv+tYyYe—-1=0, (4) 
uv + 6=0, (5) 
SO = “ ig) (6) 


(Punkt bedeutet Ableitung nach der Zeit, Strich Ableitung nach der Koordinate x, 
Index 0 Gleichstromanteil). 
Dazu die Kopplungsgleichung 


I,=i+C,E (7) 
mit 
dU,=Edx. (8) 


C, = F/(4 a) ist die Langskapazitat der Verzégerungsleitung pro Langeneinheit, 
uber welche Kreis und Elektronenstrahl miteinander gekoppelt sind. Gleichung (7) 
entspricht der fundamentalen Maxwellschen Gleichung iiber die totale divergenz- 
freie Stromdichte 


, fe 
Gee Uvaerasancn 


(In unseren Gleichungen sind Stré6me statt Stromdichten und Ladungen pro 
Lange 1 statt Ladungsdichten verwendet.) Es ist von ausschlaggebender Bedeu- 
tung, dass £ das totale Feld ist und sich in ein inneres £; und ein ausseres E, auf- 
_ spalten lasst (vgl. Arbeit [1)). 
E, rihrt von der Raumladung her und ist gegeben durch: 
. 1 
E,=— CG 
iMieeoe—— 0 Keine st elix)) 
E,, ist das Feld, das durch die Ladungen auf der Kopplungskapazitat Q = I;, 
fiir 7 = 0 (kein Elektronenstrahl) entsteht. 


. tf 
Ei = C," 


Die eingangs erwahnten Erweiterungen werden nun einzig durch den Ersatz von 


Gleichung (7) durch 
T= a OC, EB} (7’) 


das heisst im allgemeinen Fall fiir mehrere Elektronenstrahlen, erreicht. Wie in 
der Arbeit [I] werden aus den Gleichungen (1) bis (6) und (7’) einerseits das 
influenzierte Feld E in Funktion der Elektronenstréme 7, und andererseits die 
modulierten Elektronenstréme 7, als Funktion des Feldes E ausgerechnet, wobei 
E explizit in E; und E, aufgespalten ist; man erhalt bei Verwendung des Ansatzes 
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fiir alle Wechselgréssen : 


~~ et ot, (p =} @) 
yi pr a 8 
(2.-$2)-sre= Paap ey (8) 


e/M + Vox/Vor * P(r _ DZ te’ 9 
eS (p+ I tay)! (Es pC, ) } ©) 


2. Lésung der Gleichungen (8) und (9) 


Fiir die verschiedenen Fragestellungen kénnen die Gleichungen einzeln oder 
in Kombination diskutiert werden. i 
a) Die Plasmafrequenz ergibt sich als direktes Nebenresultat aus (9) fiir ein 
ruhendes Elektronenplasma v, = 0 ohne elektrischen Kreis E, = 0 mit p = 7 2 zu 


und ist schon friiher aus 4hnlichen Grundgleichungen abgeleitet worden?). 
b) Ein Elektronenstrahl. Fiir vy + 0 erhalt man aus (9) mit dem Wert (10) die 
charakteristische Gleichung fiir I’ 
Q2 
1 poe 2 SE eee (11) 
(P + I v9)? 

c) Zwei Elektronenstrahlen (Doppelstrahlverstarker). Es folgt analog aus (9) 

durch Elimination von 7, und 2, 
We tee ae Fe 2 = 

(pb + I'v)? (P + I’ vo)? 


(12) 


d) Ein Elektronenstrahl mit Helix (Travelling-Wave-Tube). Hierfiir ergibt die 
Elimination von 7 und E = E, — (i/p C,) aus (8) und (9) 


= ES een 


[a — fis (p+ I'v)? 


= Oi (13) 


I) ist der Wert fiir die Ausbreitungskonstante, den man aus (8) fiir die elektro- 
nenstrahlfreie Helix 1 = 0 erhalt (vgl. Arbeit [I)). 


3. Diskussion und Vergleich mit anderen Resultaten 


Den Vergleich machen wir mit einer von R. WARNECKE, O. DOHLER und 
W. KLEEN gemachten Publikation?), 

b und c stimmen mit a und b der Arbeit [II] tiberein. Die Raumladungskrafte 
und die Wechselwirkungskrafte zwischen zwei Strahlen erscheinen hier als Wir- 
kung des inneren Feldes gemiss Gleichung (7). 


1) G. Joos, Lehrbuch der theoretischen Physik, 5. Auflage (Akademische Verlagsgesellschatt 
Becker & Erler, Leipzig 1943), S. 676. 


*) Wireless Engineer 28, 167 (1951), Part. 1, Nr. 6 (hier zitiert: Arbeit {I1}). 


a om te gre 
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Beziiglich d tritt eine Diskrepanz mit der obgenannten Arbeit und den dort 
zitierten Autoren auf. Die Beriicksichtigung der Raumladungskrafte gibt namlich 
hier kein neues Resultat gegeniiber den Berechnungen ohne Raumladungskrafte 
(Fall c in Arbeit [II], welcher gleichbedeutend mit dem Resultat der Arbeit [1] 
ist). Fall d der Arbeit [II] ist nach dieser Berechnung hinfallig. 

Der Grund scheint uns in folgendem zu liegen: Die Berechnung der anderen 
Autoren erfolgt nicht unter Bentitzung der fundamentalen Kopplungsgleichung 
(7), sondern nur mit Hilfe des ausseren Feldes, zum Beispiel mit einer Energie- 
bilanz (Arbeit [II]). Das innere Feld wird durch eine Zusatzbetrachtung einge- 
fiihrt. Dadurch entsteht eine Unsymmetrie in den Gleichungen. Wenn zum Bei- 
spiel in Gleichung (8) links & durch £, ersetzt wird, so bekommen wir mit (9) 
auch eine von 2 abhangige charakteristische Gleichung fiir die Travelling-Wave- 
Tube; sie entspricht dem Fall d der Arbeit [II] bzw. den dort zitierten Autoren 
Cuu und Jackson, RypBecK, LAposTOLLE. Der tiefere physikalische Grund 
liegt offenbar darin, dass durch das mitlaufende Feld in der Helix die Raum- 
ladungskrafte kompensiert werden, was bei freien Elektronenstrahlen nicht der 
Fall ist, und dies scheint unabhaingig vom Kopplungsfaktor a zwischen Elektro- 
nenstrahl und Helix zu sein (vgl. Arbeit [1]; « hier = 1 gesetzt). 

Die Folgerungen aus diesen Betrachtungen sind also in Ubereinstimmung 
mit den Messungen von CUTLER und QuarTeE!’) und anderen Autoren, wonach durch 
die Uberlagerung der beiden Raumladewellen 


; Q a Q 
a ~j72(1+2), Ir, = jo (1- ) 


nach (11) langs eines freien Elektronenstrahles im Triftraum Maxima und Minima 
der Verdichtungen auftreten und die Rauschzahl eines Klystrons bei passender 
Wahl des Laufwinkels zwischen Kathode und Modulationsresonator auf ein Mini- 
mum herabgedriickt werden kann. Dies gilt auch fiir den Eintritt in eine Travel- 
ling-Wave-Tube, nicht aber fiir deren Interaktionsraum. 


Summary 


It is shown that plasma oscillations have an effect in the field-free drift space 
of the klystron and in the double-beam tube, but not in the interaction space of 
the travelling-wave tube. 


(Eingegangen: 20. Mai 1952.) 


1) C. C. CurteR und C. F. Quarter, Phys. Rev. 80, 875 (1950). 
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Frithjahrstagung der Schweizerischen Physikalischen Gesellschaft 
vom 3. Mai 1952 in Baden 


Berichte tiber angewandte Physik und Mathematik 


Die Stabilitatskriterien der Elastomechanik, von H. ZIEGLER, ETH.., Ziirich. 


Die mechanischen Systeme kénnen?) in konservative [und zwar (a) nicht- 
gyroskopische und (b) gyroskopische] sowie nichtkonservative [(c) dissipative, 
(d) zirkulatorische und (e) instationare] Systeme eingeteilt werden. Die bei Sta- 
bilitatsuntersuchungen meist verwendeten statischen Verfahren (Gleichgewichts- 
und Energiemethode) sind aber nur bei Systemen der Gruppe (a) mit dem 
kinetischen gleichwertig. Bei anderen, vor allem solchen der Gruppe (d), fiihren 
sie vielfach zu Fehlschliissen?). 


Die Knicklast des tangential gedriickten Stabes, von M. Beck, Ziirich$). 


Da die Richtung der tangentialen Druckkraft durch den Ort des Stabend- 
punktes allein nicht bestimmt wird, ist diese Belastung nichtkonservativ: es 
liegt ein zirkulatorisches System vort). Die Knicklast muss deshalb mit dem 
kinetischen Stabilitatskriterium berechnet werden’); dabei zeigt sich eindriick- 
lich das Versagen des statischen Verfahrens. 


Elektronisches Rechengerat zur Auswertung von Windkanalmessunégen, 
von H. R. VoELimy, Ziirich’). 


In der Messtechnik steht man immer wieder vor der Notwendigkeit, eine 
Vielzahl von zunachst bedeutungslosen Messgréssen in endgiiltige Messresultate 
umzurechnen. Diese Auswertearbeit benG6tigt nicht nur viel Zeit, sie ist auch eine 
standige Quelle fiir Rechenfehler. Es liegt auf der Hand, hier Rechengerate ein- 
zusetzen, welche die Auswertung rasch und fehlerfrei durchfiihren. 

Im folgenden wird ein elektronisches Rechengerat beschrieben, das zur Aus- 
wertung von 3-Komponenten-Messungen mittels Straingage-Wagen am Wind- 
kanal dient. Die drei Komponenten M, W und A miissen dabei nach folgenden 
Formeln berechnet werden: 


M _ Oy My + Cg Ny 
’ 
q q ees, 
N = Cg + Cy Mo, 
Ww T cosa + N sing f N 
ee os : a =a*+ce,N 
q q , 5 , 
7 : T = Cgt+ Cre 
A N cosa—T sing : aha ie 
q q ; 


1) Wie eine noch unveroffentlichte Untersuchung zeigt. 
2) Vgl. H. Zrecier, Ing.-Arch. 20, 49 (1952). 

8) Assistenz fiir technische Mechanik, ETH. 
4) Vgl. Referat von H. Zircier, Die Stabilitdtskriterien der Elastomechanitk. 
5) Vgl. M. Becx, ZAMP 3, 225 (1952). 

6) Contraves AG. 
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wobei die Gréssen ”,, ”, und ¢ die Messgréssen sind, «* und g Einstellgréssen 
bedeuten und die Koeffizienten c,...c, feste Werte sind. Das in Figur 1 dargestellte 
Blockschema zeigt, wie in diesem Analogiegerat obige Formeln durch Addieren 
und Subtrahieren von Spannungen und durch Multiplizieren von Spannungen 
mit Wellenstellungen realisiert werden. Mittels Potentiometern werden die Mess- 
und Einstellgréssen in ihre entsprechenden Spannungswerte verwandelt, wobei ein 
auf jedes Potentiometer umschaltbares Nachlaufvoltmeter die zahlenmassige Kon- 
trolle der eingestellten Werte erlaubt. Sobald alle Eingangsgréssen eingestellt sind, 
werden mit dem gleichen Nachlaufvoltmeter der Reihe nach die Resultatwerte «, 
M/q, W/q und A/q abgetastet und wiederum durch das Zahlwerk angezeigt. 


| 


| 
| 
— 


ee 


‘Gao | 
Lo : 


Fig. 1 Fig. 2 


In Figur 2 ist das Schema des Nachlaufvoltmeters dargestellt. Die einzustel- 
lende (oder zu messende) Spannung e wird mit der Spannung an einem genau 
linearen Potentiometer verglichen (sogenannter Helipot, Prazisionspotentiometer 
mit total 10 Umgangen). Die Differenzspannung speist tiber einen Verstarker 
den Motor M stets so, dass die Differenzspannung verschwindet (sogenannter 
Nachlauf). Ein an der Potentiometerwelle mit entsprechendem Ubersetzungs- 
verhaltnis angeschlossener Umlaufzahler zeigt somit direkt das Verhaltnis der 
Spannung e zur Spannung U am Potentiometer an. Der auf der gleichen Welle 
sitzende Generator G bewirkt mit seiner Gegenspannung lediglich eine Stabilisie- 
rung des ganzen Nachlaufsystems, wobei die Nachlaufzeitkonstante am Potentio- 
meter P eingestellt werden kann. 

Der grosse Vorteil des Vergleichsvoltmeters mit Nachlauf besteht neben der 
zahlenmissigen Angabe der Werte in seiner Unempfindlichkeit gegen Schwankun- 
gen der Referenzspannung U, da stets ein Spannungsverhaltnis bestimmt wird. 

Das Rechengerat, das nicht grosser als ein normaler Radio ist, enthalt demnach 
neben dem Netzteil einen Réhrengenerator mit Verstarkern zur Erzeugung der 
Referenzspannung U und zur Speisung der Motor-Generator-Gruppe, das Nach- 
laufsystem mit den dazugehérenden Verstarkern und das Rechennetzwerk. 

Die Genauigkeit des Rechengerates hangt im wesentlichen von der Linearitat 
des Helipots ab und betragt rund 2 Promille, was fiir die meisten technischen An- 
wendungen geniigen diirfte. 

Wahrend der Messung am Windkanal konnen nun von einem Gehilfen die 
laufenden Auswertungen vorgenommen und sogleich graphisch aufgetragen 
werden, so dass der Versuchsingenieur sein Messprogramm stets den erhaltenen 


Resultaten anpassen kann! 


396 . Varia — Miscellaneous — Divers ZAMP 


In weiterer Vereinfachung des Mess- und Auswertevorganges ist vorgesehen, 
die von den Straingage-Wagen gelieferten Spannungen direkt in das Rechengerat 
zu leiten, so dass der Vorgang des Ubertragens der Werte wegfallt. Ferner werden 
die Resultate von mehreren Vergleichsvoltmetern auf Schreibgerate tibertragen, 
so dass die gesuchten aerodynamischen Beiwerte direkt als Funktion des ebenfalls 
aufgezeichneten Anstellwinkels aufgetragen werden. 


Analogierechengerate mit linearen Potentiometern, von A. P. SPEISER, 
Ziirich), 

Unter Verwendung von linear gewickelten Prazisionspotentiometern der kauf- 
lichen Bauart und von ohmschen Widerstanden lassen sich Rechengerate bauen. 
Die elementaren Rechenoperationen werden durch einfache Schaltungen ausge- 
fiihrt. Ausserdem kénnen vorgegebene Funktionen mit einer oder mehreren Va- 
riablen durch Approximation mittels rationaler Funktionen dargestellt werden; 
entsprechende Verfahren zur Ermittlung der Schaltungsparameter sind entwik- 
kelt worden. Auf dieser Grundlage sind Gerate zur Auswertung von standig wie- 
derkehrenden mathematischen Ausdriicken sowie zur Lésung von nichtlinearen 
Gleichungssystemen gebaut worden. Die Gerate kénnen mit bescheidenen Kosten 
und innert kurzer Frist gebaut werden; sie sind robust, einfach zu bedienen und 
benétigen nur einfachste Stromquellen. Ein ausfiihrlicher Bericht iiber diese 
Konstruktionen wird demnachst erscheinen. 


Schaufelschwingungen, von Cu. BELLENoT, Baden?). 


Der Autor fiihrte einen Film iiber Schwingungen an Fliigelgittern vor und 
gab dazu einige theoretische Erlauterungen. 


Eine neue Auswertemethode fiir Hitzdrahtmessungen, von G. DATWYLER, 
Ziirich. 


An Hand der Kingschen Gleichung (auch in stark verallgemeinerter Form) 
liber das Warmegleichgewicht bei Hitzdrahtsonden in Strémungen lasst sich 
zeigen, dass ein einfacher linearer Zusammenhang besteht zwischen aquivalenten 
Geschwindigkeits- und Heizstromanderungen, die entgegengesetzt gleiche Tem- 
peratur- und Widerstandsaénderungen der Sonde bewirken. 

Der Proportionalitatsfaktor zwischen den aquivalenten Geschwindigkeits- 
und Heizstromanderungen ist eine eindeutige und universell giltige, ein fiir alle- 
mal berechenbare Funktion der Reynolds-Zahl des Hitzdrahtes in der Strémung, 
wie durch zahlreiche Messresultate experimentell bestatigt ist. 

Hieraus lasst sich eine neue, einfache und rasche Auswertemethode fiir Tur- 
bulenzmessungen mit dem Hitzdrahtanemometer gewinnen, 


Wahl der Grundbahnen zur Bestimmung ganzer Flugbahnscharen bei 
Raketen, von E. Rorn, Luzern, und R. SANGER, ETH.., Ziirich. 

In einer friiher erschienenen kurzen Mitteilung*) wurde darauf hingewiesen, 
wie ein Satz von PoINncAR# iiber die Parameterabhangigkeit der Lésungen von 
Differentialgleichungen zur Berechnung von Flugbahnscharen herangezogen 
werden kann. Sollen bei der Interpolation die extremalen Abweichungen von den 
wahren Werten dem Betrage nach einander gleich sein, so bedeutet dies eine spe- 


1) Institut fiir Angewandte Mathematik der ETH. 
2) AG. Brown, Boveri & Cie. 


3) E. Rotu-DESMEULEs, Zur Berechnung der Flugbahnscharen ferngesteuerter Raketen, ZAMP 2, 
487-489 (1951). 
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zielle Wahl der Parameterwerte der zu berechnenden Grundbahnen. Es wird 
gezeigt, dass nach geeigneter Variablentransformation diese Werte sich als Null- 
stellen von Tschebyschewschen Polynomen ergeben. 


Synchronisierung von Reflexoszillatoren, von A.H. ABDEL DAYEM, Ziirich!). 


Es ist bekannt, dass man einen selbsterregten Generator durch eine einge- 
pragte Fremdspannung synchronisieren kann, sofern die Frequenz des Fremd- 
signals angenadhert in einem einfachen rationalen Verhaltnis zur Frequenz des 
selbsterregten Signals steht. Es existiert immer eine bestimmte Synchronisie- 
tungsbreite, welche sowohl von der Grésse des Fremdsignals und den Daten des 
Schwingsystems als auch von den Eigenschaften der nichtlinearen Kennlinie des 
Generators abhangig ist. 

Wie in den verschiedenen Frequenzgebieten der Hochfrequenztechnik findet 
man auch im Mikrowellengebiet wichtige Anwendungen fiir den synchronisierten 
Generator. Daher ist die Kenntnis von Theorie und Anwendung der Synchroni- 
sierung selbsterregter Generatoren auch im Mikrowellengebiet von grosser Wich- 
tigkeit. In einer demnachst erscheinenden Arbeit®) wird diese Frage insbesondere 
am Reflexoszillator (Reflexklystron) untersucht, derim Bereiche von 8 cm bis 1 cm 
Wellenlange Verwendung findet, wenn keine grossen Leistungen erforderlich sind, 
wie zum Beispiel beim Uberlagerungsempfang, bei Signalgeneratoren fiir Mess- 
zwecke usw. ; 

Im wesentlichen besteht ein Reflexoszillator aus einem Schwingsystem (Hohl- 
raumresonator), das durch einen gesteuerten Strom erregt wird, das heisst einen 
Strom, der von der Spannung abhangig ist, die am Schwingsystem liegt. Diese 
Strom-Spannungs-Abhangigkeit ist durch eine nichtlineare Funktion darstellbar, 
die durch den Mechanismus von Geschwindigkeitsmodulation und .«Bunching» 
gegeben ist. Mit Hilfe dieses vereinfachten Bildes lasst sich der eingeschwungene 


2S) 


Pigzel 
Verhalten des durch ein Fremdsignal synchronisierten Reflexoszillators. 


Iz Ausgangsleistung im synchronisierten Zustand; 
P, Ausgangsleistung im «ungestérten» Zustand; 

ora optimale Ausgangsleistung vom Reflexoszillator ; 
P, Leistung des eingepragten Fremdsignals; 

0) Verstimmung} 

| to) | 7 Grenzverstimmung des Synchronisierungsgebietes. 


1) Institut fiir Hochfrequenztechnik der ETH. 
2) Dissertation ETH. (1952). 
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Zustand des synchronisierten Reflexoszillators berechnen. Fir einen durch ein 
unabhangiges Fremdsignal synchronisierten Reflexoszillator ergibt die theoretische 
Berechnung die Resultate in Figur 1. 

Fir den Fall der gegenseitigen Synchronisierung zweier identischer Reflex- 
oszillatoren sind die theoretisch berechneten Kurven in Figur 2 dargestellt. 


Fig. 2 
Verhalten zweier identischer Reflexoszillatoren im gegenseitigen Synchronisieren. 


P  Ausgangsleistung von einem Reflexoszillator; 

Pi, Leistung gekoppelt von einem Oszillator zum anderen; 
Oy, relative Verstimmung der Hohlraumresonatoren ; 

wy Lange des Kopplungsweges, dargestellt als Winkel. 


Das Experiment bestatigt in allen Teilen die Aussagen der Theorie. 


Ein Laboratoriumsgerat zum Auffinden von Undichtigkeiten in Hoch- 
vakuumsystemen (Leak-detector), von G. WEIBEL, Ziirich), 


Angesichts der technologischen Schwierigkeiten, hochvakuumdichte Gefasse 
und Systeme (Glas-Metall-Verschmelzungen, Schweiss- und Létnahte, Gussteile 
usw.) herzustellen, haben Gerate zur Lokalisierung von Undichtigkeiten grosse 
Bedeutung gewonnen. Die Empfindlichkeit solcher Instrumente ist je nach Me- 
thode und apparativem Aufwand sehr unterschiedlich. 

Im Institut fiir Hochfrequenztechnik an der ETH. wurde fiir das Arbeiten 
mit demontierbaren, experimentellen Réhrensystemen und verschiedenen Hilfs- 
apparaten das nachfolgend kurz beschriebene Gerit entwickelt, welches hohe 
Empfindlichkeit bei relativ geringem Aufwand erzielt. Das Priifen findet unter 
normalen Vakuum-Betriebsbedingungen statt. Mit einem Suchstrahl von Testgas 
wird das System von aussen her abgestrichen, waihrend ein dem Vakuumraum 
angeschlossener Indikator auf Spuren eindringenden Testgases reagiert. 

Das System spricht auf den beim Uberstreichen der Undichtigkeit hineindif- 
fundierenden Testgasimpuls genau so an wie ein elektrisches Netzwerk aus 
R- und C-Elementen auf einen eingespiesenen Stromstoss (Analogien in Grund- 
gleichungen) ; Zeitkonstanten, Impulsflankensteilheit und asymptotischer Druck- 
anstieg lassen sich aus den Systemgréssen ermitteln. 

Die Verhaltnisse an der Indikatormesszelle, hier beim vorliegenden Gerat ein 
normales Magnetron-Ionisationsmanometer («Philips-Manometer»), wurden ge- 


1) Institut fiir Hochfrequenztechnik der ETH. 
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nauer betrachtet ; sowohl die Anderungen der Partialdriicke Luft—Testgas als auch 
damit verbundene Anderungen des Totaldruckes bewirken eine differentielle 
Variation des Manometerstromes. In einer Briickenschaltung oder mit Hilfe eines 
differenzierenden Netzwerkes wird sodann eine Signalspannung gewonnen, einem 
streng symmetrischen Gleichstromverstarker und Modulator zugefiihrt, welcher 
ein akustisches Signal mit einer der Eingangsspannung proportionalen Frequenz 
erzeugt. 

Theoretisch kann der gebaute Leak-detector bei einem Vakuumrezipienten 
von V= 51 und einer effektiven Saugleistung S = 2 1/s beispielsweise im Druck- 
gebiet von p = 10-* bis 10-* mm Hg eine Einrinnmenge von OQ 1OO WS marin Lele 
mit einer Zeitkonstanten von T= 2,5 s noch deutlich anzeigen (df = 100 Hz). 
Durch interne Gasquellen (Gasausbriiche diverser Art bei nicht gut entgasten 
Systemen) wird beim beschriebenen Indikator die nutzbare Empfindlichkeit herab- 
gesetzt. Es wurde auf eine symmetrische Briickenschaltung mit zwei Manometer- 
réhren verwiesen, die es gestattet, mit einfachen Mitteln zwischen Testgas und 
systemeigenem Gas zu diskriminieren (gasselektive Anzeige) und damit die volle 
Empfindlichkeit auszunutzen. 

Literatur: 

F. M. Pennine, Tiefdvuckmanometer, Philips’ techn. Rdsch. 2, 201-208 (1937). 
F. M. Pennine und K. Nirenuuis, Bauart und Anwendung einer neuen Ausfih- 
vung des Philips-Vakuummeters, Philip’s techn. Rdsch. 77, 116-123 (1949). 


tiber den Einfluss der Plasmaschwingung in Laufzeitrohren, von ip, Lupi, 
Baden). 


Es wird gezeigt, dass die Plasmaschwingung im feldfreien Triftraum des 
Klystrons und in der Zweielektronenstrahlréhre einen Einfluss hat, nicht aber im 
Interaktionsraum der Travelling-Wave-Tube. Ausfiihrliche Publikation siehe 
ZAMP 3, 390 (1952). 


Warmespannungen in elastischen Kérpern, von F. Sarzmany, Ziirich?). 
Ausfiihrliche Darstellung dieses Problems siche ZAMP 3, 129 (1952). 


Der Einfluss der Molzahlainderung auf den Wirkungsgrad von Ver- 
brennungsmaschinen, von W. TRAUPEL, Winterthur?). 


Es wird gezeigt, dass die Anderung der Molzahl, die bei der Verbrennung im 
allgemeinen auftritt, eimen Einfluss auf den Wirkungsgrad von Verbrennungs- 
maschinen (Verbrennungsmotoren und Gasturbinen) ausiibt. Molzahlvergrosse- 
rung verbessert den thermischen Wirkungsgrad, Molzahlverkleinerung verringert 
ihn. Am Beispiel eines idealisierten Gasturbinenprozesses wird dies anschaulich 
gemacht. Bei Molzahlvergrosserung ware es grundsatzlich sogar denkbar, aus 
einem Kreisprozess Arbeit zu gewinnen, selbst wenn die ganze fiihlbare Warme 
(Heizwert) dem Prozess sogleich wieder nach aussen entzogen wird. Damit besteht 
auch grundsatzlich die Moglichkeit, den Carnot-Wirkungsgrad zu iiberschreiten. 
Fs wird dabei eine gewisse Warmemenge der ausseren Umgebung entzogen und 
in Arbeit verwandelt, was keinen Widerspruch zum zweiten Hauptsatz darstellt, 
da ja ausserdem eine chemische Umsetzung vollzogen wird und die Entropie- 
summe im ganzen zunimmt. 


1) AG. Brown, Boveri & Cie. 
2) Forschungsabteilung Escher-Wyss AG. 
3) Gebr. Sulzer AG. 
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Bei technischen Prozessen wird der Einfluss in bestimmten Fallen deutlich 
bemerkbar, und zwar im negativen Sinne (Molzahlverminderung bei Brenngasen, 
die einen grossen Gehalt an CO aufweisen). Ein Brennstoff, bei dem ein grosser 
positiver Effekt erzielt werden kénnte, ware Methylalkohol. Genaueres siehe: 
W. TRAUPEL, Der Einfluss des Brennstoffes auf den Wirkungsgrad von Verbren- 
nungsmaschinen, Allg. Warmetechn. 3, 1-9 (1952). 


Rechnungs- und Messergebnisse iiber Hautkondensation an gewellten 
Oberflachen, von R. Grecoric, Ziirich?). 

Nach einer einleitenden Ubersicht iiber Haut- und Tropfenkondensation wird 
die Méglichkeit der Hautkondensation an einer unebenen Oberflache erértert. 
Die Oberflachenspannungen rufen, der verschiedenen Kriimmungen der Konden- 
satoberflache wegen, sehr starke Druckgradienten hervor. Letztere kénnen ein 
grosses Vielfaches (ein paar Zehnerpotenzen) grésser als diejenigen des Schwere- 
feldes sein. Starke Druckgradienten ergeben, aus Gleichgewichtsgriinden, diinne 
Kondensathaute, was nach der Theorie von NuSsELT zu grossen Warmeiiber- 
gangszahlen fiihrt. Fiir einen Kondensator mit vertikalen Rohren, die mit Rillen 
parallel zur Rohrachse versehen waren, sind fiir gegebene geometrische Verhalt- 
nisse die mit Differenzenrechnung erhaltenen Warmetibergangszahlen angege- 
ben worden. Letztere sind durch Versuche bestatigt. Zwei Kennzahlen der Modell- 
ahnlichkeit bestimmen die Nusseltsche Zahl des Warmeiibergangs. An Hand 
dieses Ahnlichkeitsgesetzes, das mit Versuchen bestatigt worden ist, wurden fiir 
verschiedene Dampfe die Warmeiibergangszahlen an gerillten Rohren mit jenen 
an glatten Rohren verglichen. Dabei betrug das Verhaltnis der letztgenannten 
Warmeiibergangszahlen etwa (5-7) zu 1. 


Paul Felix Neményi 
June 5, 1895 — March 1, 1952 


NemEny1’s life-long interest was pure mechanics. He worked mainly in five 
fields: static elasticity, fluid dynamics, hydrology and hydraulics, organization 
of mechanics, methods of research. ‘ 

The first of his three major discoveries in elasticity was a new singularity 
method [Z. angew. Math. Mech. 9, 488 (1929); 70, 383 (1930)]. A number of rather 
complicated singularities, obtained by confluence of simple ones, are classified 
in terms of their influences in accordance with a duality principle, which gene- 
ralizes MAXWELL’s reciprocity theorem. 

His continuing interest in analogies may well have begun with his own 
discovery that the stream-lines of any potential flow of an incompressible fluid 
may serve also as stress trajectories for a plane elastic system. He first approached 
the problem through his method of singularities (Proc. int. Congr. appl. Mech. 
1930); later [Z. angew. Math. Mech. 72, 364 (1931)], he gave an analytic proof 
of the beautiful reformulation now called “NEMENYI’s Theorem”: given any 
net of isothermal curves, there exists a five parameter family of plane stress systems 
for which these curves are stress trajectories, Phrased thus in terms of pure elasticity, 
it suggests another question: to what extent is the solution of a problem in con- 
tinuum mechanics characterized by an associated trajectory system? This sub- 
ject attracted Nem&ny1 throughout the rest of his life. 


1) Escher-Wyss, Maschinenfabriken AG. 


Vol. III, 1952 Buchbesprechungen — Book Reviews — Notices bibliographiques 401 


His third major discovery in elasticity is the reduction of the general exten- 
sional theory of thin shells of revolution with meridian y = /(¥) to the remarkably 
simple equation 


From ‘“‘ NEMENyt1’s stress functions” U,, the general solution is obtained imme- 
diately [Proc. 6th int. Congr. Math. (1936), Proc. nat. Acad. Sci. 29, 159-162 
(1943)]. 

All this theoretical work illustrates the inverse or semi-inverse approach, 
which he summarized in a fine organizational summary [Advances in Applied 
Mechanics, Vol. II (Academic Press, New York, 1951), p. 123-151]. In papers 
written jointly with Prim [Proc. nat. Acad. Sci. 34, 119 (1948); J. Math. Phys. 
27, 130 (1948) and others], he applied it to obtain numerous exact solutions of 
the nonlinear equations of gas dynamics, many of them representing rotational 
and three-dimensional flows of nonuniform total energy. 

For any person with a genuine question in mechanics, NEMENYI was always 
willing to supplement his great knowledge and deep understanding by hours 
of library work. Since 1946 many of his ideas were worked out by his pupils (in 
temporal order): C. TRUESDELL, R. C. Prim, A. VAN Tuvt, A. W. SAENZ, R. Tou- 
PIN, J. L.EricKsen. He had held posts in many institutions in Europe and 
America; at the time of his death he was Head of the Theoretical Mechanics 
Section at the U. S. Naval Research Laboratory. 

C. TRUESDELL, Indiana University. 
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Die Fakultat (Gammafunktion) und verwandte Funktionen, mit be- 
sonderer Beriicksichtigung ihrer Anwendungen. Von L6scH-SCHOBLIK 
(B. G. Teubner, Verlagsgesellschaft, Leipzig 1951). 205 S., 22 Abb.; $4.03. 


Die von L. EuLer in die Analysis eingefiihrte Gammafunktion (Fakultat) 
und die verwandten Funktionen kommen in so vielen Gebieten der Mathematik 
vor und sind fiir so zahlreiche Anwendungen ein unentbehrliches Instrument, 
dass die Kenntnis ihrer Haupteigenschaften zum Riistzeug des Mathematikers 
und des Physikers gehdren muss. Das vorliegende Buch — dessen Bearbeitung 
durch F. Léscu unter Zugrundelegung eines von F. ScHoprik hinterlassenen 
Manuskriptes erfolgte — gibt eine Einfithrung in die Theorie und in die Anwen- 
dungen dieser Funktionen. Die Theorie ist aber an Formeln und Entwicklungen 
so reichhaltig, dass eine Auswahl nétig war; sie wurde im Hinblick auf die An- 
wendungen auf Physik und Technik getroffen. Das Studium des Werkes kann 
jedem empfohlen werden, der im Besitze der Elemente der Theorie der analyti- 
schen Funktionen ist und sich eine griindliche Kenntnis der Gamma- und der 
verwandten Funktionen aneignen will. 

Inhaltsangabe: I. (S. 1-100). — 1. Die Produktdarstellungen und die funda- 
mentalen Eigenschaften der Fakultat. — 2. Die asymptotischen Darstellungen der 
Fakultat (Stirlingsche Formel). — 3. Die Reihendarstellungen der Fakultat. — 
4. Die Integraldarstellungen der Fakultat. — 5. Auswertung bestimmter Integrale 
durch Fakultaten. — 6. Die Funktionalgleichung R,(z) f(z + 1) + R2(2) f(z) = 0. 


ZAMP III/26 
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Anwendungen auf die hypergeometrische Reihe und auf das Integral von MELLIN- 
BARNES. 


II. (S. 101-151). — 1. (z, g)! und Q(z, @) als Funktionen von z. — 2. (2, o)! und © 


Q(z, g) als Funktionen von g. — 3. Die Integralexponentiale und verwandte Funk- 
tionen (Integralkosinus und -sinus usw.). — 4. Das Fehlerintegral und verwandte 
Funktionen. — 5. Die Hermiteschen Polynome. — 6. Die Funktionen des para- 
bolischen Zylinders und die Whittakerschen Funktionen. 

III. (S. 152-203). Anwendungen. — 1. Eine elektrostatische Aufgabe. — 2. Die 
Sitzreiheniiberhéhung im Theater. — 3. Zur strahlentheoretischen Treffertheorie. 
— 4, Uber kontinuierliche Leibrenten. — 5. Der Strahlungswiderstand einer Ver- 
tikalantenne. — 6. Wahrscheinlichkeits- und Fehlerrechnung. — 7. Warmeleitung. 
— 8. Die astronomische Strahlenbrechung. — 9. Beugung elektromagnetischer 


a 


Wellen. — 10. Ubergangsbégen bei Bahnen. M. Plancherel — 


The Lebesgue Integral. By J.C. BurxKitt, Cambridge Tracts in Mathe- 
matics and Mathematical Physics, No. 40 (Cambridge University Press, Cam- 
bridge 1951). 87 pp.; 12s, 6d. 


Depuis la parution en 1905 du livre de LEBESGUE Legons sur l’intégration et 
la vecherche des fonctions primitives de nombreuses et importantes simplifications 
ont été apportées par DE LA VALLEE PoussIN et d’autres a l’exposition et a la 
démonstration des théoremes. I] est maintenant possible d’exposer cette théorie 
sous une forme extrémement simple et accessible a toute personne tant soit peu 
familiarisée avec les éléments du calcul différentiel et intégral. Celle que présente 
M. BurkILt dans cet opuscule est basée sur la théorie de la mesure des ensembles; 
elle répond pleinement aux conditions de clarté et de simplicité que l’on peut 
maintenant exiger et sa lecture peut étre recommandée a qui désire étre mis 
rapidement au courant des théoremes essentiels de la théorie de l’intégrale. 

Contenu: I. Ensembles de points. II. Mesure. III. L’intégrale de LEBESGUE. 
IV. Différentiation et intégration. V. Autres propriétés de l’intégrale. VI. L’inté- 
grale de LEBESGUE-STIELTJES. Solution des exemples et des problémes. 

M. Plancherel 


Einfihrung in die theoretische Physik. Bd. 3, Teil 2, Quantentheorie. 
Von C. SCHAEFER (W. de Gruyter, Berlin 1951). 510 S., 88 Abb.; DM 40.—. 


Der in der zweiten Auflage vorliegende Band III, 2 enthalt die Lehren der 
Quantentheorie bis und mit den Diracschen Gleichungen. Zunachst werden die 
Schwierigkeiten der klassischen Physik in der Theorie der Strahlung, der spezi- 
fischen Warme, des Photoeffekts usw. geschildert. Es folgt an Hand der Bohr- 
schen Theorie eine sich tiber vier Kapitel erstreckende Darlegung des Atom- 
modells von Bour und die sich daraus ergebenden Folgerungen fiir Atombau, 
Spektren und periodisches System. Dann wird in Kapitel 6 die Wellenmechanik 
eingefiithrt im Anschluss an DE BRoGLie und ScHrOpDINGER. Ihre Anwendung 
und Deutung behandeln Kapitel 7 und 8. Darauf wird nochmals das Problem 
der Strahlung aufgegriffen und, jetzt mit wellenmechanischen Methoden, ein- 
gehend geschildert (Auswahl- und Polarisationsregeln, Ubergangswahrscheinlich- 
keiten, Dispersion und Ramaneffekt). Das zehnte Kapitel schliesslich enthalt die 
Diracsche Theorie des Elektrons und Positrons. 

Gegeniiber der ersten Auflage sind keine wesentlichen Anderungen im Sinne 
einer Modernisierung vorgenommen worden, was im Hinblick auf die erprobte 
Darstellungsweise des Verfassers zu bedauern ist. W. Baumgartner 
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Funzioni ellittiche. Di FRANcEsco Tricomi (N. Zanichelli, Bologna 1951). 
353 p.; lit. 4500.—. 


L’auteur s’est proposé de donner au lecteur qui connait les premiers éléments 
de la théorie des fonctions analytiques un exposé aussi simple que possible de la 
théorie des fonctions elliptiques lui permettant de l’appliquer sans difficulté a des 
problemes concrets. Laissant de cdté les applications a l’algébre et a la théorie 
des nombres, auteur a réussi a rassembler dans cet ouvrage l’essentiel de ce que 
doit connaitre celui qui veut appliquer les fonctions elliptiques aux problémes 
posés par la physique et la technique et pousser leur solution jusqu’aux résultats 
numériques. L’ouvrage répond donc pleinement au but cherché. 

Un premier chapitre expose la théorie des fonctions elliptiques de WEIER- 
STRASS et la représentation d’une fonction elliptique générale par ces fonctions. 
Il se termine par l'étude du comportement des fonctions de WEIERSTRASS a in- 
variants réels. 

Le second chapitre traite des intégrales elliptiques, de leur réduction a trois 
espéces fondamentales et aux formes canoniques de LEGENDRE. I] contient un 
tableau trés pratique des formules qui ramenent une intégrale elliptique aux 
formes de LEGENDRE dans le cas ot le polynome du 3¢ ou du 4°¢ degré a ses 
coefficients réels. 

Le chapitre III introduit les fonctions elliptiques de JacoBr et les fonctions 
théta; il les relie aux fonctions de WEIERSTRASS et montre l’utilité des fonctions 
théta pour le calcul numérique des intégrales elliptiques de 3¢ espéce. 

Le chapitre IV est consacré a l'étude des transformations du 1° et du 2° ordre 
des fonctions de WEIERSTRASS et de JACOBI. 

Le dernier chapitre, d’environ 80 pages, expose lusage des tables numé- 
riques des fonctions et des intégrales elliptiques et donne de nombreuses appli- 
cations 4 des problémes tels que la rectification de l’ellipse et de hyperbole, la 
recherche des géodésiques de l’ellipsoide de rotation, le flambement d’une poutre, 
le pendule simple, la détermination ‘d’une transformation conforme a propos d’un. 
probléme aérodynamique. II est suivi d’une table récapitulative des principales 
formules établies dans le cours de l’ouvrage (18 pages), d’un tableau comparatif 
des notations de l’auteur et de celles employées dans les ouvrages classiques, 
d’une liste bibliographique et d’un indice alphabétique. 

A noter pour le lecteur auquel la langue italienne n’est pas familiere qu’une 
traduction allemande de la premiére édition par M. KrarFr a paru en 1948 


a Leipzig (Akad. Verlagsgesellschaft) sous le titre Elliptische Funktionen. 
M. Plancherel 


Lectures on Classical Differential Geometry. By D. J. StrurK (Addison- 
Wesley Press, Inc., Cambridge 1950). 221 pp., 120 figs.; $6.00. 

Das vorliegende Werk ist aus einer einsemestrigen Vorlesung des Verfassers 
am Massachusetts Institute of Technology hervorgegangen und legt die Grund- 
begriffe der klassischen Kurven- und Flachentheorie dar. Die einzelnen Kapitel 
sind iiberschrieben mit: Kurven, Elementare Flachentheorie, Die Fundamental- 
gleichungen der Flachentheorie, Geometrie auf der Flache, Spezielle Gegenstande 
aus der Kurven- und Flachengeometrie (Enveloppen, Abbildungen und Land- 
karten, Minimalflachen, Das Imaginare in der Flichentheorie). 

Das Buch ist sehr ausfiihrlich gehalten. SrrurK verwendet die Vektorschreib- 
weise, wie sie heute fiir diesen Gegenstand fast allgemein iiblich ist. Die notwen- 
digen Hilfsmittel aus der Vektorrechnung sowie aus andern Zweigen der Mathe- 


matik sind jeweils, bevor sie gebraucht werden, in Kleinschrift in den Text ein- 
: @ 
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gebaut. Die einzelnen Kapitel sind durch Aufgaben abgerundet, von denen die 
meisten vom Leser ohne fremde Hilfe gelést werden kénnen. So kann sich der 
Leser iiberpriifen; gleichzeitig offenbaren aber diese Aufgaben 6fters einen inter- 
essanten geometrischen Sachverhalt. Eine Serie von weniger elementaren Auf- 
gaben ist am Schlusse des Buches zusammengestellt. 

Besondere Aufmerksamkeit hat der Verfasser den Illustrationen geschenkt; 
sie zeichnen sich durchgehend durch grosse Anschaulichkeit aus. Einige davon 
sind dusserst gut geratene Photographien von mathematischen Modellen. Ebenso 
verdienen die zahlreichen historischen und biographischen Notizen Beachtung. 

Das Buch iiber Differentialgeometrie von STRUIK kann jedermann empfohlen 
werden, der sich iiber diesen Gegenstand orientieren méchte. Es eignet sich be- 
sonders zum Selbststudium. M. Jeger 


Integraltafel: Erster Teil: Unbestimmte Integrale; Zweiter Teil: Bestimmte 
Integrale. Herausgegeben von W. GROBNER und N. HOFREITER (Springer-Verlag, 
Wien und Innsbruck, 1949/50). VIII + 166 S. bzw. VI + 204 S. sFr. 23.50 und 
sFr. 25.-. 


Die beiden Tafeln enthalten zu den angegebenen Integralen auch theoretische 
Anleitungen fiir den Gebrauch; ferner sind in vielen Fallen Rekursionsformeln 
angegeben, die entlastend wirken. Bemerkenswert ist die gegeniiber ahnlichen 
Tafeln etwas starkere Betonung der Bediirfnisse der reinen Mathematik. 

1. Teil: Unbestimmte Integrale. Beginnend mit den Integralen von rationalen 
Funktionen mit Erlauterungen zur Partialbruchzerlegung bietet das Buch weiter- 
hin eine reichhaltige Sammlung von Integralen, die Wurzeln enthalten. Insbe- 
sondere sind Anleitungen zur Umrechnung von elliptischen Integralen auf die 
Legendresche Normalform gegeben. Schliesslich folgen noch transzendente Inte- 
granden in mannigfacher Form und am Schluss des ersten Bandes Integrale iiber 
Jacobische und Weierstrasssche elliptische Funktionen. 

2. Teil: Bestimmte Integrale. Dieser zweite Band enthalt eine Fiille von wert- 
vollem Material, insbesondere zahlreiche Integrale, die nur zwischen bestimmten, 
numerisch gegebenen Grenzen in geschlossener Form ausgewertet werden kénnen. 
Der Praktiker wird eine solche Tafel um so mehr schatzen, als bestimmte Integrale 
in der Literatur meist schwer zu finden sind. Im iibrigen sind die verschiedenen 
Typen von Integralen nach den gleichen Gesichtspunkten angeordnet wie im 
ersten Band, jedoch mit einem besonderen Abschnitt iiber bestimmte Integrale 
von Zylinderfunktionen. 

Die besondere Reichhaltigkeit dieses Tabellenwerks macht es zu einem wert- 
vollen Hilfsmittel fiir das Integrieren. H. Rutishauser 


Proceedings of a Second Symposium on Large-Scale Digital Com- 
puting Machinery (Harvard University Press, Cambridge, USA., 1951). 293 pp., 
121 figs., $8.-. 

Das Werk (Band 26 der «Annals of the Computation Laboratory of Harvard 
University») ist ein ausfiihrliches Protokoll eines Kongresses an der Harvard 
University im September 1949, Die 38 gehaltenen Vortrage sind wiedergegeben. 
Die mathematischen wie die technischen Aspekte programmegesteuerter Rechen- 
maschinen sind in gleicher Weise beriicksichtigt. 

Unter den technischen Arbeiten finden wir Beschreibungen einiger im Betrieb 
oder im Bau befindlicher Rechenautomaten, ferner Erlauterungen von neuen 
elektrischen Schaltelementen; Erwahnung verdienen hier die Aufsatze Static 
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Magnetic Delay Lines, Coordinate Tubes for Use with Electrostatic Storage Tubes, — 
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Electrochemical Computing Elements. Weitere Beitrage betreffen neue mathema- 
tische Methoden; insbesondere sind folgende Arbeiten hervorzuheben: Numerical 
Methods Associated with Laplace’s Equation, An Iteration Method for the Solution 
of the Eigenvalue Problem of Linear Differential and Integral Operators, On the 
Monte Carlo Method. Andere Arbeiten behandeln sodann Anwendungen von pro- 
grammgesteuerten Rechenmaschinen auf Probleme der theoretischen Physik, 
der Aerodynamik, der Statistik und der Nationalékonomie. Den Abschluss bildet 
ein bemerkenswertes Referat The Future of Computing Machinery von L. N. 
RIDENOUR. 

Das Buch stellt eine interessante Dokumentation tiber den Stand des pro- 
grammgesteuerten Rechnens im Jahre 1949 dar und vermittelt dem Fachmann 
zahlreiche wichtige Anregungen. Ae S CUS CV 


Grenzschicht-Theorie. Von H. ScHLicHT1ING (G. Braun, Karlsruhe 1951). 
483 S., 295 Abb.; DM 45.-. 


Die Strémungslehre erhielt einen entscheidenden Impuls, als PRANDTL 1904 
seine Grenzschichttheorie aufstellte. Zum ersten Male war es méglich, die bei 
grésseren Reynoldsschen Zahlen auftretenden Widerstande auf Korper zu ver- 
stehen. Dadurch, dass die Reibung vor allem in einer relativ diinnen Oberflachen- 
schicht wirkt, wahrend die 4ussere Str6mung praktisch reibungsfrei lauft, wird 
die ungeheuer schwierige Aufgabe der gleichzeitigen exakten Beriicksichtigung 
von Reibung und Tragheit in zwei mindestens naherungsweise viel leichter lés- 
bare Teilaufgaben zerlegt. Schon friih gelang es Brasrus, das Plattenproblem 
exakt zu lésen. Ein weiterer wichtiger Schritt gelang 1913/14, als man den Unter- 
schied von laminarer und turbulenter Grenzschicht bei reibungs- und turbulenz- 
freier Aussenstrémung klar erkannte. Von da an ist ein stetiger Fortschritt zu 
verzeichnen, nachdem dann ausser der Géttinger Schule sich auch die wbrige 
Welt an der Grenzschichtforschung beteiligte. 

Die Grenzschichttheorie wird in allen Strémungslehrbiichern behandelt. Es 
fehlte aber bisher (mindestens im deutschen Sprachgebiet) an einer Monographie, 
die alles Wesentliche iiber Grenzschichten zusammenfasst und damit das mitih- 
same Suchen in der verstreuten Literatur erpart. SCHLICHTING, der selbst funda- 
mental wichtige Beitrage vor allem zum Problem der Stabilitat laminarer Schich- 
ten geleistet hat, war wie kaum einer berufen, diese Aufgabe zu iibernehmen. So 
liegt denn ein stattlicher Band vor, dessen Studium jedem, der naher mit Grenz- 
schichten zu tun hat, aufs warmste empfohlen werden kann. 

Zuerst wird die Zahigkeit in elementarer Weise eingefiihrt und ihre Wirkung 
beschrieben. Nachher folgen die exakten Gleichungen und diejenigen Falle, die 
exakt gelést werden kénnen. Entsprechend dem Ziel des Buches sind die Falle 
mit starker Reibung (Stokes-Formel, Lagerschmierung) nur kurz behandelt. Es 
folgt dann die ausfiihrliche Behandlung der laminaren Grenzschichten — ohne und 
mit Druckgradienten der Hauptstromung; auch sind nichtstationare Falle be- 
handelt. Die Naherungsmethoden sind naturgemass sehr eingehend dargestellt, 
wobei sowobl ebene als auch rotationssymmetrische Falle gleichmassig behandelt 
werden. Besonders interessant ist das Kapitel iiber Grenzschichtbeeinflussung, 
wo ziemlich viele neue Ergebnisse, die in der Kriegs- und Nachkriegszeit gewonnen 
wurden, mitgeteilt werden. Kompressible laminare Grenzschichten mit und ohne 
Warmeiibergang folgen. 

Das nun folgende Kapitel iiber die Entstehung der Turbulenz scheint dem 
Referenten eher etwas knapp gefasst. Man mochte gerne, dass der Verfasser, der 
ja gerade auf diesem Gebiet zu Hause ist, die Methodik der Stabilitatsrechnungen 
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etwas ausfiihrlicher darstellen wiirde. Hingegen ist die Ubersicht iiber die bis- 
herigen Lésungen sehr ausfiihrlich und nitzlich. 

Der voll turbulenten Grenzschicht und den freien Strahlen ist der ganze Rest ; 
des Buches gewidmet. Auch hier ist die neueste Entwicklung eher etwas zu kurz 
gekommen; der mehr elementare Standpunkt, der dem Ingenieur natiirlich be- 
sonders niitzlich ist, wiegt hier vor. Besonders interessant ist auch das Kapitel 
iiber Rauhigkeit. 

Man muss das Schlichtingsche Buch als ein Standardwerk bezeichnen, das 
fiir jeden, der mit der Strémung wirklicher Fliissigkeiten und Gase zu tun hat, 
eine unentbehrliche Fundgrube darstellt. Ganz besonders werden die Ingenieure 
dariiber gliicklich sein, da es der Verfasser verstanden hat, von den mathemati- 
schen Mitteln weisen Gebrauch zu machen, und die physikalische Wirklichkeit 
nicht aus den Augen verloren wird. J. Ackeret 


Nomography and Empirical Equations. By L. Jonnson (John Wiley 
& Sons, New York, 1952). 150 pp.; $3.75. 


Dieses neue Lehrbuch behandelt zwei verschiedene, verwandte Gebiete. Der 
erste, 90 Seiten umfassende Teil ist der Nomographie gewidmet, der zweite, 
kiirzere, der analytischen Approximation empirischer Funktionen. 

Im nomographischen Teil erértert der Verfasser zunachst eingehend den 
Begriff der Funktionsleiter. Sodann entwickelt er die Konstruktionsregeln der 
Fluchtentafeln mit drei parallelen Leitern, der Z- oder N-Tafeln und der Tafeln 
mit zwei parallelen und einer krummlinigen Leiter, welche zu den Funktions- 
typen f, + fe = fs, f1 = fo’ fs und f, - g3 + fe + fg = 0 gehGren. Diese Tafeln ent- 
halten bekanntlich Funktionsleitern, die in f/, und f, regular sind. In gewissen 
Fallen, etwa wenn diese Funktionen sehr grosse Werte annehmen, kann es wiin- 
schenswert sein, in bezug auf 1/f/, und 1/f, regulare Leitern zeichnen zu kénnen. 
Dies kann durch geeignete kollineare Abbildung der erwahnten Tafeln erreicht 
werden. Es ist ein Verdienst des Verfassers, diese Abarten von Tafeln und Funk- 
tionstypen besonders behandelt zu haben. — Eingehend werden auch Beziehungen 
zwischen vier und mehr Veranderlichen erfasst, die durch Verbindung der bisher 
genannten Tafeltypen mittels Zapfenlinien darstellbar sind. 

Mit diesen einfachen Mitteln beherrscht der Autor eine grosse Menge prakti- 
scher Beispiele aus der Technik, der er auch zahlreiche Aufgaben entnimmt. Auf 
kompliziertere Tafeltypen verzichtet er, und er lasst, was wesentlicher ist, auch 
das ganze Gebiet der Netztafeln beiseite. Insofern ist der Titel des Buches zu 
weit gefasst. 

Im zweiten Teil (Empirical Equations) geht es darum, beobachtete Folgen 
von Wertepaaren einer Funktion zu approximieren mittels einer Funktion y (x) 
von vorgeschriebenem analytischem Bau, in der Hauptsache Potenzfunktionen 
oder weniggliedrige Polynome und Exponentialfunktionen. Als Verfahren zur 
Berechnung der Konstanten wird zwar die Methode der kleinsten Quadrate 
skizziert. In beinahe allen Beispielen wird jedoch untersucht, ob das Beobach- 
tungsmaterial bei Eintragung auf geeignetes Funktionspapier (einfach- oder dop- 
peltlogarithmisch usw.) sich geradlinig anordnet, die Funktion also hier ersten 
Grades wird. Durch Riickkehr zu den Variablen x und y ergibt sich dann das 
Gesetz y(*). 

Der Schreibende méchte an JoHNSoNS Werk die leichtfassliche und elemen- 
tare Darstellung loben, besonders aber das klare Herausarbeiten der Resultate 
und ihre Zusammenstellung in Tabellen, die dem Beniitzer sehr niitzlich sein 
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Die Messwandler - ihre Theorie und Praxis 


Von Dr. J. Gotpstern, dipl. Ing. 


Das Buch stellt eine vollstandig neu bearbeitete und wesentlich erganzte zweite Auflage 

des im Jahre 1928 unter dem gleichen Titel erschienenen Werkes dar. Der Entwick- 

lung auf dem Messwandlergebiet und den wichtigsten Neuerscheinungen der letzten 
Jahre wird Rechnung getragen. 


222 Seiten mit 199 Abb. In Ganzleinenband Fr. 29.10, broschiert Fr. 24.95 (1952) 


Die Gestalt der elektrischen Freileitung 


Von Dr. techn. MILAN VIDMAR 


Der vom Transformatorenbau her sehr bekannte Verfasser, der sich seit ungefahr 

einem Jahrzehnt der elektrischen Energietibertragung zugewandt hat, versucht in 

diesem Buch die Grundziige einer Theorie aufzubauen, die die bisherige ziemlich 
rohe, obgleich erfolgreiche Empirie des Freileitungsbaues ersetzen soll. 


200 Seiten mit 49 Figuren. Ganzleinen Fr. 19.75, broschiert Fr. 16.65 (1952) 
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Anwendung e Wirkungsweise e Schutzmassnahmen 
U. JETTER unter Mitarbeit von W. HANLE 


Unter Verwendung des offiziellen Berichts der U.S. Atomic Energy Com- 

mission versucht diese Verdffentlichung, von wissenschaftlicher Seite her 

das Grundsitzliche der Wirkung von Atombomben und die Méglichkeiten 

von Schutzmassnahmen darzustellen. Sie soll dazu dienen, der bewusst gross- 

gezogenen Angst vor der Atombombe entgegenzutreten und eine richtige 
Einschatzung dieser Waffe zu gewinnen. 


Aus dem Inhalt: Die Atombombe als Objekt der Politik. — Physikalische und 
chemische Grundlagen: Kernprozesse; Uran und Plutonium; Kettenreaktion 
im Kernreaktor und in der Atombombe; Messverfahren der Kernphysik. - 
Eigenschaften und Einsatz von Atomwaffen: Atombomben und radioaktive 
Gifte; Die wirkliche und die nominelle Atombombe; Die Luftexplosion ; 
Explosion im Wasser, auf dem Boden und in der Erde. —: Die Wirkung von 
Atombomben: Druckwelle und Luftstoss; Licht- und Warmestrahlung; Kern- 
strahlung und radioaktive Verseuchung. — Planung von Schutzmassnahmen : 
Die Erfahrungen in Japan; Allgemeine Schutzmassnahmen ; Kernstrahlung 
und radioaktive Verseuchung; Planung auf weite Sicht. 


72 Seiten A5, 30 Abbildungen und Zeichnungen, 8 Tabellen, brosch. DM 4.80 
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J.J. BURCKHARDT 


Professor an der Universitat Ziirich 
p - 


Die Bewegungsgruppen der Kristallographie | 


nf 3 
(1947) 184 Seiten mit 56 Figuren. Leinenband Fr. 33.30, broschiert Fr. 29.10 q 


| 


«Das Ziel des Buches ist eine rein mathematische und in allen Einzelheiten auch 4 
dem Nichtmathematiker verstandliche Behandlung der Kristallklassen und Bewe- 
gungsgruppen. Die friihere einschlagige Literatur war meistens entweder nur dem _ | 
Mineralogen oder nur dem Mathematiker zuganglich. Es gelang dem Verfasser, | 
den Anspriichen beider Arten von Lesern zu entsprechen, indem er alle mathema- 
tischen Hilfsm'ttel und Uberlegungen klar darlegt, nur durch die kristallographische 
Zielsetzung bedingte mathematische Hilfsmittel gebraucht und die (Schoenflies- q 
schen) Bezeichnungen und Benennungen der Kristallographie beniitzt, andrerseits _ 
aber seinen Stoff in rein mathematischer Weise behandelt und so auch das Inter- 
esse des reinen Mathematikers fiir den aus der Kristallographie entsprungenen | 
Gegenstand erweckt. » (Zentralblatt fiir Mathematik) — | 
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VERLAG BIRKHAUSER - BASEL UND STUTTGART 


Neuerscheinungen aus dem Verlag Birkhauser, Basel und Stuttgart 


Lehrbticher und Monographien aus dem Gebiete der exakten Wissenschaften 
Mathematische Reithe — Band 13 und 18 


L. BIEBERBACH 


Theorie der geometrischen Konstruktionen 
(1952) 162 Seiten mit 102 Figuren. In Ganzleinenband Fr. 18.70, broschiert Fr. 15.60 


Das Buch ist aus den Vorlesungen des Verfassers an den Universitaten Basel, Frankfurt 
am Main und Berlin entstanden. Es gibt eine systematische Darstellung von Problemen 
und Ergebnissen, die seit dem Altertum in Forschung und Lehre ein stetes Interesse gefunden 
haben. Die durch die neuere Literatur und durch eigene, hier erstmals veréffentlichte Uber- 
legungen des Verfassers méglich gewordenen Fortschritte wurden eingearbeitet. 


H. HERRMANN 


Ubungen zur projektiven Geometrie 


(1952) 168 Seiten mit 90 Figuren, Betrachtungsbrillen und 4 Raumbildern 
In Ganzleinenband Fr.17.—, broschiert Fr. 14.— 


In den vorwiegend analytischen Lésungen zu 279 Ubungen wird auch der Kenner einige 

neue Darstellungsweisen finden, so in der Art der Verwendung der Matrizen, deren for- 

male Eleganz durch inhaltlich-anschauliche Deutungen bereichert erscheint. Die Ubungen 

sind zum Gebrauch neben Vorlesungen und Lehrbiichern bestimmt, insbesondere ergan- 
zen sie die Biicher von W. BLAscuKE in dieser Sammlung. 
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- Mathematische Werke aus dem Verlag Birkhdnser Basel / Stu 


Der Briefwechsel von Johann Bernoulli. Herausgegeben von der Bernoulli-Kommission der 
_ Naturforschenden Gesellschaft Basel. Band I: ca. 480 Seiten mit zahlreichen Figuren . . . 


_ _Brepersacu, L.: Theorie der geometrischen K onstruktionen. 170 Seiten mit 103 Figuren . . . 


_ Buascuke, W.: Analytische Geometrie. Zweite, verbesserte Auflage. ..... - . ate 
_ Biascuxe, W.: Projektive Geometrie. Dritte, verbesserte Auflage. . - . . . Pee et 
_ Burcxuarpt, J. J.: Die Bewegungsgruppen der Kristallographie ............ 
- Caratutopory, C.: Funktionentheorie. In zwei Banden. 


Band I: 288 Seiten mit 33 Figuren. - . ...... +s. .-+5 5. AO OR hs 
Batid 11; 194 Seiten mit 73. Figuren j°-04 7. 6 see cee oe ee eee See Sa ke be 


. ‘Finsier, P.: Uber Kurven und Flachen in allgemeinen Raumen. Unveranderter Nachdruck 
b _ der Dissertation von 1918, vermehrt um ein ausfiihrliches Literaturverzeichnis von 
Dr Hs Schubert: 170,Seiten ie) Sy 2°.) A/S (ae cae eden ee AD we Le 


Hurwitz, A.: Mathematische Werke 
_ Band 1: Funktionentheorie. 734 Seiten mit 238 Figuren. . .........2.2.4... 


Liner, A.: Statistische Methoden fiir Naturwissenschafter, Mediziner und I) ngenteure. 150 Seiten 
mit 38,Figuren. Zweite, erweiterte Auflage .) )o@w a) .)6 1. to eon eee 


Locuer-Ernst, L.: Differential- und Integralrechnung im Hinblick auf ihre Anwendungen. y To 
596 Seiten mit 406 Figuren. . . . . . 5 fe Te ete ee io se 3G 8 ed a re i (isa wae eels 49.90 


Locuer-Ernst, L.: Einfihrung in die freie Geometrie ebener Kurven. 85 Seiten mit 168 Fig. brosch. 12.50 
Ostrowskl, A.: Vorlesungen tiber Differential- und Integralrechnung 4 ’ yes 
_ BandI: Funktionen einer Variablen. 373 Seiten mit 42 Figuren. . ...... ‘ geb. 56.1! 


_ Band II: Differentialrechnung auf dem Gebiete mehrerer Variablen. 484 Seiten, 55 Fig. geb. 69.70 
Band III in Vorbereitung. ari 


” 


TEFEL, E.: Lehrbuch der darstellenden Geometrie. 174 Seiten mit tiber 180 Figuren . . ... eb. 29.65 


VorxeR, D., und Doetscu, G.: Die sweidimensionale Laplace-Transjormation. 260 Seiten — “yo 
mit 17 Figuren und vielen Tabellen ©. 2°. -s)ee ty 6 sp Chee ee geb. 47.85 
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Archiv der Mathematik — Archives of Mathematics — Archives des mathématiques oo 

_ Herausgegeben in Verbindung mit dem Mathematischen Forschungsinstitut in Oberwolfach von 
H. Kyeser und W. Svss. : 

Beirat: G. Bor, E. Bompian1, P. Ten Bruccencate, J. Diruponn&é, Cu. EHRESMANN, H. GORTLER, 

H. Hapwicer, H. Horr, W. Manus, T. NaGELL, Cur. Pauc, J. Rapon, K. REIDEMEISTER, 

. J.A. Scuouren, H. Srirert, R. Sperner, E. STiEFeL. Redaktion: H. Bituarz, Universitat Freiburgi. Br. 


Erscheint seit 1948; ab Band 8 (1952) jahrlich in 6 Heften von je ca, 80 Seiten Umfang. Abonnementspreis 
pro Band Fr. 60.—, Einzelheft Fr, 12.—. ; 
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_ Zeitschrift fiir angewandte Mathematik und Physik (ZAMP) ; a 
Journal of Applied Mathematics and Physics — Journal de Mathématiques et de Physique appliqueés — 
Redaktionskommission: J. Ackrrrt, E. BAUMANN, P. Nicci, P. SCHERRER, E. St1EFEL, F. Srissi, _ 


H. ZinciEr. Redaktion: R. SANGER, Eidg. Technische Hochschule, Ziirich. Erscheint seit 1950 alle 2 
4 Monate. 16 X 23 cm. Jahrlich Fr. 34.—; Ausland Fr. 40.—; Einzelnummer Fr. 8.—. 
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